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ABSTRACT  (ONGERUBRICEERD) 


The  GLUSE  model  has  been  developed  to  describe  die  dynamics  of  (fresh  water)  qpray  droplets  in  die 
hmnogeneous  boundary  layer  of  an  air-sea  interaction  simulation  tunnel,  and  dieir  influence  on  humidity  and 
tenqierature  profiles  piouault  et  aL,  1991].  ExtensitHi  to  opoi  ocean  conditions  requires  that  the  fiosh  water 
droplets  are  rqilaced  by  sea-salt  aerosol,  and  diat  die  influence  of  waves  is  taken  into  account.  In  this  rqiort 
we  describe  the  progress  on  the  development  of  the  SeaQuse  model  that  t^lies  to  open  ocean  conditions, 
llie  calculation  of  die  profiles  of  die  mean  wind  velocity  and  turbulent  diifusivity  has  been  changed  to  take 
into  account  the  charactoistics  of  the  marine  atmosphere  and  of  die  wave  field.  A  module  was  added  to  the 
initialization  routines  diat  models  die  air  flow  over  the  waves  in  a  non-turbulent  and  non-eviqKnative 
atmoqihere.  This  module  yields  vertical  profiles  of  droplet  concentrations  that  enter  in  the  relaxatkm  term 
diat  models  the  mean  movonent  of  die  particles.  These  profiles  are  dien  modified  in  the  iterative  budget 
calculation  diat  accounts  for  droplet  tranqxHt  atmosidietic  turbulence  and  dieir  interaction  with  die 
humidity  field  by  evaprvation.  In  this  part  of  the  code,  die  evqKmtkm  module  was  changed  to  describe  die 
evrqiaratioo  of  sea  salt  aerosols  in  die  marine  atmo^diere.  The  ouqxit  of  die  new  SeaQuse  model  consists  of 
vertical  profiles  of  tempoature,  humidi^  and  aerosol  concentration  from  the  sea  surfiroe  iqi  to  the  top  of  the 
marine  atmoqiheric  boundary  layer. 
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SAMENVATTING  (ONGERUBRICEERD) 

Het  CLUSE  model  werd  ontwikkeld  om  de  dynamica  van  (zoet  water)  qjray  diuppel^  in  de 
1-dimensionale  bomogene  grenslaag  van  een  hidit-water  interactie  tunnel  le  beschrijven,  alamede  de  invloed 
van  de  waterdtuppel^  op  de  vodttig|>eids-  en  tmqteratuuipiofielen  [Rouault  et  al.,  1991].  Voor  de 
uidneiding  van  het  model  naar  'open  ocean'  condities  is  vereist  dat  de  zoet  water  diuppels  vervangen  worden 
door  zeezout  aerosolen,  en  dat  de  invloed  van  de  golven  woidt  gemodelleerd.  In  dit  lappcnt  wordt  de 
voortgang  besclireven  van  het  geschikt  maken  van  h^  CLUSE  model  voor  toepassing  in  'open  ocean* 
condities.  De  berekening  van  de  verticale  profielen  van  de  gemiddelde  windsnelbeid  en  de  turtnilente 
difiiisie  coefficknt  zijn  aangqiast  aan  een  maritienie  atmosfeer  en  golven.  De  inirioliMtfe  van  het  modd  is 
nitgebieid  met  een  module  waarin  de  hichtstroining  over  de  golven  wordt  besclireven  in  een  atmosfeer, 
waarin  geen  tuibolentie  of  verdamping  voodDomt  Deze  module  berdcent  de  concentratie  profiden  vm  het 
aerosol  die  gdimikt  worden  in  de  relaxatie  term  die  de  gemiddelde  beweging  van  de  deel^  besduijft.  Deze 
ctmcentratie  profiden  wnden  vervolgens  gewijzigd  in  een  iteratieve  budget  berekening,  waarbij  tranqxHt 
van  het  aerosol  door  atmosferisdie  tuibolentie  en  verdaii^ping  van  het  aerosol  wordt  gemodelleerd.  De 
module  die  de  venlanq;ring  van  de  aerosol  deel^  bescfarijft,  werd  aangepast  voor  zeezout  aerosolen  in  de 
maiiene  atmosfeer.  De  resoltaten  van  bet  nieuwe  SeaChue  mcdd  bestaan  uit  de  verticale  profielen  van 
tenywatour,  vochtigheid  en  aerosol  concentratie  vanaf  het  zee  oppetvlak  tot  aan  de  top  van  de  maiiene 
grenslaag. 
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1  INTR(M>UCn<»J 

Aerosols  play  an  important  role  in  die  marine  boondaiy  layer,  bodi  in  meteorological  processes 
and  in  die  field  of  propagatum  of  electro-optical  radiation.  Suspended  aerosols,  1^  scattering  and 
absorption  of  radiation,  contribute  to  die  attenuatUMi  of  radiation  in  die  electromagnetic 
transmission  windows.  This  reduces  die  performance  of  electro-optical  systems  and  limits  electio- 
optical  surveillance  of  die  marine  boundary  li^er  and  die  sea  surfsce  from  aircrafts  and  satellites. 
The  reduced  albedo  may  also  be  of  inyortance  for  climatological  studies. 

On  a  micro-metenological  scale,  aerosols  contribute  to  die  moisture  and  heat  budgets  near  die  sea 
surface  [Wu,  1974;  Ling  and  Kao,  1976;  Andreas,  1992].  Since  die  ocean  acts  bodi  as  a  source 
and  a  sink  for  atmoyberic  aerosols,  die  aerosol  dnylets  may  transfer  water  vapor,  heat,  poUutants 
and  bacteria  diroogh  the  air-sea  interfroe.  At  higher  levels,  sea-salt  particles  shrink  by  evaporation 
or  grow  by  condensation  trough  interaction  with  die  humidity  field.  Hence,  diey  play  a  role  in  die 
transport  of  matter,  heat  and  water  vqior  duough  the  marine  boundary  layer  and  are  important  to 
larger  scale  meteorological  processes  and  climatology.  Li  additkm,  small  droplets  may  act  as 
condensation  nuclei  and  are  dietefore  inycatant  in  die  formation  of  fog  and  clouds. 

In  view  of  the  above,  adequate  descr^ons  of  die  particle  size  distribotions  for  the  ambioit 
atmospheric  conditions  ate  required  to  eiylain  a  wide  range  of  processes  in  the  marine 
atmos{dieric  boundary  layer.  The  Navy  Aerosol  Model  (NAM)  [Gathman,  1983]  is  (me  of  die 
most  widely  accqpted  and  used  models  to  assess  die  propagatkm  diaracteristics  of  electro- 
magnetic  radiatkm.  It  is  die  maritime  versicm  of  die  LOWTRAN  ctxle  [Knei^  etJd.,  1983]. 
NAM  is  an  enyirical  modeL  based  <m  an  extensive  data  set  including  a  variety  of  geognqdiical 
locations.  Over  die  years,  NAM  has  been  lydated  and  improved  widi  bodi  eiqierimental  [De 
Leenw,  1986]  and  dieoretical  [Gerber,  1985]  evklence  [Hughes,  1987;  Gathman,  1989]. 

NAM  lyplies  to  a  single  hei^t  ctf  10  m  above  mean  sea  level  (AMSL).  Therefore,  aikiitional 
modeling  is  retpured  to  infer  die  vertical  profiles  of  aerosol  ctmcentration.  To  diis  end  die  Navy 
Oceanic  Vertical  Aerosol  Modd  [Gathman,  1989;  De  Leenw  etal.,  1989a]  is  being  devdcyed.  hi 
NOVAM,  die  concentratkm  of  aerosol  at  deck  height  is  provided  by  NAM,  and  die  vertical 
distributioo  of  aerosols  is  subsequently  calculated  widi  physical  equations. 
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Althoag^  NAM  and  NOVAM  provide  na  widi  valuable  mfonnatinn,  they  cannot  be  uaed  below 
10  m  AMSL.  Yet,  diia  regum  ia  important,  becanae  an  appieciaUe  part  of  die  pnqMgatkm  pad>  for 
detecdon  of  targeta  low  above  the  horizon  be  below  10  m  AMSL.  The  region  ia  alao 
in^ortant  for  climatological  atudiea,  becanae  the  cmqiling  tenna  between  atmoqdiere  and  ocean 
in  global  climate  modela  are  not  yet  ascertained.  As  described  above,  an  inqxMtant  part  of  the 
aerosol  dynamics  alao  takes  place  in  the  sniface  layer.  Therefore,  any  model  that  is  developed  for 
dm  marine  atnio8|dieric  boundary  layer  requites  a  detailed  desaqrtkm  of  die  processes  near  die 
air-sea  interface. 


The  processes  near  the  air-sea  interface  are  difGcnlt  to  assess,  because  diey  consist  of  a  large 
number  of  coiqiled  dynamical  and  j^ysical  processes.  When  wave  breaking  occurs,  air  is 
entrained  into  the  water.  The  air  breaks  iq;>  in  bubbles,  which  rise  due  to  dieir  buoyancy.  When 
they  reach  the  surface  and  burst,  sea-spray  aerosol  is  produced  (see  figure  1.1).  Sea  quay  conaiats 
of  small  film  droplets,  which  result  from  the  lueaking  of  die  diin  water  film  covering  the  bobble 
when  it  protrudes  die  water  surface,  and  jet  druplets,  which  result  from  the  btealdng  op  of  die  jet 
that  is  sobseqoentty  formed. 


■  —  ■■  ■—  ..  ,  .,^z=sh 

Turbulent 


Figniel.l:  of  breaking  waves  with  bobUe  and  spray  drofriet  formation  (from  [LarsmetaL, 
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The  production  of  film  and  jet  droplets  from  boA  single  bubbles  and  bubble  plumes  has  been  the 
subject  of  extensive  ejqierimental  roidies  (see  e.g.,  Blandtard  [1963,  1983],  Mtmahan  etal. 
[1982],  Nfesti^er  and  Lefimconnier  [1988],  Spiel  [1992a]).  An  analytical  descripticn  of  die  jet 
droplet  production  has  recently  been  provided  [Dekker  ami  De  Leeuw,  1993],  baaed  on  analysis 
by  Machityre  [1968, 1972].  In  addidmi  to  die  film  and  jet  droplets,  the  so-called  qnime  droplets 
are  jsoduced  by  wave  crest  tearing  udien  die  wind  spttA  exceeds  9  m/s  [Monahan  etaL,  1986; 
Wu,  1990;  Andreas,  1992]. 


The  number  and  size  distribution  of  sea  qny  aerosols  diat  is  produced  is  given  by  die  droplet 
source  function.  Monahan  et  aL  [1982]  presented  a  source  function  diat  consists  of  two  modes, 
diat  is,  a  bubble-mediated  and  a  qnime  droplet  source  function.  Alteniative  formulations  have 
been  given  by  Fairall  etal.  [1983]  (see  also  Miller  and  Fairall  [1988]),  De  Leeuw  [1990a]  and 
Andreas  [1992]. 


WIND  - ► 


MIXING  LAYER 


Sdieniaticrqwesentationofdie  wave-rotor  model 
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However,  it  is  very  difBciilt  to  obtaiii  pertinent  measmements  of  die  iqiwanl  suifKe  fhn  of 
droplets  at  sea  and  die  extension  of  laboratory  lesnlts  to  oceanic  conditions  is  not  stndg^itfotward 
either.  As  a  result,  diere  is  not  yrt  a  consensus  source  function  availabk.  In  particular,  die 
contzibotion  of  spume  droplets  has  by  no  means  been  established  [Katsaros  and  De  Leeuw,  1993]. 

The  fieshly  produced  droplets  are  tranqiorted  iqiward  by  turbulent  diffusion  and  convective  flow, 
counterbalanced  by  gravitational  forces.  While  they  are  aiibotne,  diey  intenux  widi  the  scalar 
fields  of  tenqierature  and  humidity  by  beat  exchange,  ev^xnatum  and  condensational  growdi. 

Fbw  measuremrtits  exist  of  the  concentration  of  aerosols  below  ship  deck  level  (see  De  Leeuw 
[1989b,  1993]  for  a  review).  Until  recently  die  particle  concentratioiis  were  assumed  to  increase 
exponentially  toward  the  sea  surface  [Blanchard  and  Woodcock,  1980].  De  Leeuw  [1986,  1987] 
provided  e^qieiimental  evidence  that  die  profiles  are  not  always  logaridumc  and  hypothesized  a 
wave-rotor  model  to  eiqilain  diis  behavior  (see  figure  12).  According  to  diat  modeL  the  aeosols 
are  tenqiorarily  tnqqied  in  eddies,  «bich  results  in  longer  suspension  times.  Alternatively,  Wu 
[1990]  suggested  diat  die  spume  droplet  production  by  die  wave  tearing  mechanism  was  to  be 
held  responsible  for  die  non-logaridimic  profiles. 

This  discussion  shows  dutt  diere  is  a  requirement  for  a  model  diat  gives  a  ibysical  description  of 
die  behavior  of  marine  spray  droplets  in  the  atmospheric  surfiu«  layer  close  to  die  waves,  and  diat 
provides  an  adequate  explanation  for  the  (qwrse)  eiqierimental  observations.  Ova*  the  years, 
various  interesting  models  have  been  reported  (see  Rouault  et.al.  [1991]  for  an  overview),  but  a 
detailed  descr^on  of  die  behaviour  of  droplets  close  to  die  surface  has  not  yet  been  given. 

fin  1981,  the  Humidity  Exchange  Over  die  Sea  (HEXOS)  program  was  initiated  to  study  die  water 
vipor  and  droplet  fluxes  firom  sea  to  air  (see  Smith  etal.  [1990]  for  a  descr^on).  The  HEXOS 
studies  in  a  Simulation  Tunnel  (HEXIST)  subprogram  was  aimed  at  die  development  of 
numerical  models  to  study  die  turbulent  transport  of  evqwtating  drcplets  and  dieir  interactions 
widi  the  scalar  fiekls  of  teaperature  and  humidity.  The  measurements  n^cfa  would  provide  die 
basis  for  die  model  fbrmulatkms  were  performed  in  die  period  1985-1990  in  the  Large  Air-Sea 
hiteraction  Simulation  Tunnel  of  die  histitnt  de  Mdeanique  Statistique  de  la  TVnbulence  (IMST) 
laboratory  at  Luminy  (MarseUle,  Rrance),  and  in  die  ubiternp  simuiation  tank  at  die  University  of 
Connecticut  (see,  e.g.,  Mestayer  etaL  [1988]). 
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It  was  recognized  diat  breaking  waves  m  a  laboratory  tunnel  cannot  provide  a  good  simulation  of 
die  net  production  of  aerosols  aid  water  vapor  over  ocean  waves  [Mestayer  and  Lefimcotmirr, 
1988].  Therefore,  die  tunnel  was  tqiaated  in  an  artificial  configuration,  where  die  surface  sources 
of  droplets  were  provided  by  bubble  jdumes  generated  by  inunersed  aeradon  devices,  where  there 
were  no  waves  and  where  die  air  flow  was  a  well-known  fully  turbuloit  boundary  layer  over  a  flat 
water  surface. 


Two  different  types  of  models  were  developed: 

(a)  A  Lagrangian  Monte  Carlo  ^pe  simuladon  of  droplet  trajectories  [Edson,  1989]; 

(b)  An  Eulerian  K-diSiisivity  model  [Rouault  etal.,  1991]. 

Althou^  fundamentally  different,  given  die  same  initial  condidons,  the  results  fiom  the  two 
models  show  excellent  agreement  between  each  other  and  with  die  tunnel  measurements.  We  will 
not  further  discuss  Edson's  model  in  diis  report,  but  restrict  ourselves  to  die  Eulerian  model  (b). 

The  Couche  Limite  Unidimensionelle  Stadonnaiie  d'Embcuns  (CLUSE)  model  by  Rouault  et  aL 
considers  fresh  water  jet  droplets  between  5  and  105  pm  in  radius,  ejected  fixim  die  water  surface 
by  bursting  bubbles.  The  model  is  tuned  by  an  eiqierimeatal  source  fiincdon.  The  droplets  are 
then  diqiersed  in  a  fully  developed  boundary  layer  turbulence  over  a  flat  air-water  interface. 
The  difiusion  by  air  turbuleiice  is  rqaesented  by  a  K-diffnsivity  term  and  gravitadonal  and 
merdal  effects  are  explicidy  modeled.  The  drqilet  population  is  qilit  in  size  bins  (categories). 
Each  droplet  category  is  allowed  to  interact  widi  die  ambient  humid  air.  Source-sink  fiinctums  in 
die  budget  eqputtions  fm  die  mass  concentration  of  eadi  droplet  categcny,  for  die  water  vi^ 
concentratum,  and  for  the  sensible  beat  rqxesent  die  effects  of  droplet  evqxxation:  transfer 
between  droplet  categories  due  to  shrinking,  production  of  water  vqxu  and  absorption  of  sensible, 
heat  A  brief  summary  of  die  original  CLUSE  hypotheses  and  model  equations  is  presented  in 
chapter  2  of  this  rqiort 

The  C1..USE  model  could  successfully  demonstrate  sane  of  die  mechanisms  of  die  non-linear 
intoactkm  between  injection,  evqxxation,  turbulent  difiusion  and  tranqiott  of  ^my  drpidets 
[Mestiqrer,  1990;  Rouault  etaL,  1991].  However,  since  it  was  developed  for  artificial  coulitions 
in  a  wind-wave  interaction  tunnel  (see  dxxve),  the  CLUSE  model  most  be  extended  to  open  ocean 
ctuiditions  before  it  can  describe  the  bdunriw  of  marine  qny  dnqilets. 
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The  fiirdier  devek^unent  of  die  CLUSE  model  for  qipljcatioa  over  sea  is  a  joint  ^ort  of  the  TNO 
Pl^sics  and  Electrooka  Laboiatofy  and  the  Ecole  Ccntrale  de  Nantes.  Initially,  the  model  code 
was  optimized  and  mwte  mofe  flexible  to  allow  introdoction  of  open  ocem  omditions  [De 
Leeuw,  etal.,  1991].  The  exteosim  of  CLUSE  to  SeaChise  then  requites  three  basic  steps.  One 
has  to  provide: 

(a)  a  model  for  the  air  flow  over  the  waves 

(b)  an  evqioration  module  for  sea  water  droplets 

(c)  a  source  function  for  open  ocean  omditions 

The  present  efforts  are  mainly  aimed  at  (a)  and  (b).  We  describe  die  air  flow  over  die  waves  a 
model  that  leads  to  die  wave-rotors  proposed  by  De  Leeuw  [1986].  The  eviqxiration  of  salt  water 
will  be  treated  in  die  way  that  was  outlined  previously  by  Mestayer  [1990].  Salt  water  evqxnation 
differs  from  fresh  water,  because  die  droplets  can  only  shrink  to  a  minimum  size  determined  by 
their  salt  content  and  ambient  relative  humidity. 

Chsqiter  3  presents  our  model  for  the  mean  air  flow  over  the  waves,  which  yields  vertical  profiles 
of  U(z),  W(z),  py(z)  and  T(z)  in  die  non-evaporative  and  non-turbulent  atmosj^iere.  The 
movement  of  (sea  ^nay)  droplets  in  die  non-turbulent  air  flow  and  in  die  absence  of  eviqioraticm 
is  discussed  in  cluqiter  4.  fri  die  budget  equations  of  die  SeaQuse  model  the  terms  diat  teprcaeat 
the  fluxes  of  droplets  due  to  their  average  movement  are  modelled  by  their  vertical  profiles  diat 
would  be  observed  in  the  absence  of  turbulence  and  of  evsporatimi. 

During  the  actual  run  of  die  SeaCluse  modeL  die  atmospheie  is  turbulent  and  ev^xnaticm  of 
droplets  takes  place.  The  algoridun  for  the  ev^ioration  of  die  salt  water  droplets  is  outlined  in 
chiqxer  S,  vdiereas  chiqiter  6  discusses  the  inqilemoitation  of  die  evqioration  module  in  the  Sea¬ 
Cluse  code. 

The  results  are  brkfly  discussed  in  cluqiter  7,  where  we  also  indicate  the  line  of  wok  to  further 
develop  and  test  die  SeaCluse  modeL 
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2  THE  ORIGINAL  CLUSE  MODEL 

The  CLUSE  model  was  developed  by  Rouault  and  Mestayer  and  has  been  extensively  described 
[Rouault,  1989],  [Rouault  and  Larsen,  1990],  [Rouault  etaL,  1991].  Therefore,  (xily  a 
compilatitHi  of  the  main  hypotheses  and  model  equations  is  given  here. 

The  framework  of  the  model  is  defined  by  the  following  hypodieses: 

1.  The  fluid  is  incompressible 

2.  The  Coriolis  force  is  negligible 

3.  There  are  no  chemical  reactions 

4.  The  Boussinesq  iq>proximaticHt  may  be  applied 

5.  The  droplets  are  spherical  and  do  not  split 

6.  The  number  of  droplets  in  die  domain  of  study  is  ccmserved 

7.  The  boundary  layer  is  fully  developed  and  tibe  classical  horizrmtal  honxigeneity  hypothesis  is 
iqiplied.  The  constant  flux  hypothesis  can  be  used  in  all  the  domain  for  the  transportable 
variables  that  have  no  local  source 

8.  The  density  of  droplets  is  too  small  to  have  any  influence  on  the  dynamics  of  air 

9.  Only  jet  droplets  are  taken  into  account 

The  fluid  is  considered  to  be  a  multiphase  mixture  of  N-i-2  conqionents:  dry  air,  with 
concentration  p,,  water  vrqior  with  concoitration  p„  and  N  droplet  categories  of  radius  r^  (n  > 
1.  and  concentration  Pg.  Category  n  includes  all  droplets  of  radius  r  defined  by:l‘l 

In  -  Arg/2  <  r  s  Tg  ♦  Ar^/2 

The  instantaneous  budget  equation  of  mass  of  any  of  die  N-l-2  transportable  conqionents  is: 

-  3  fT-a,v.n 

®Cj  ’  I  n  -  1, . .  ,N 


1 


Sx  appenlU  A  for  i  UMiiig  of  all  tyidbols  aaed  la  dSi  r^oit 
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wfaeie  is  Ae  veloci^  cooqKxient  in  direcdon  j  of  the  component  y  and  Sy  is  the  net  source 
tenn.  In  die  CLUSE  model,  the  exchange  of  mass  between  the  con^onents  can  only  be  due  to 
dn^let  evqxjradon.  Therefore: 


( 


S.  -  0 


(2.3) 


Hie  sl4>  velocity  is  introduced  as  die  velocity  difference  between  a  component  and  the  local 
center  of  mass.  We  use  the  decomposition  of  a  quantity  O,  which  characterizes  the  turbulent  flow, 
into  a  mean  part  4  and  a  fluctuating  part  such  that  O'  ^  0.  With  this,  and  hypotheses  7  and  8, 
the  instantaneous  budget  equatimis  can  be  avenged  to  yield  die  mean  budget  equatimis; 


3K 

“5r 


a(pE) 

dz 

3 


dz 

4  5 


(2.4) 


where  W^  and  are  the  vertical  conqioDents  of  droplet  mean  and  slip  velocity,  respectively,  and 
w'  denotes  the  fluctuating  conqwnent  of  die  vertical  air  flow.  In  arklition,  die  budget  equation  of 
sensible  heat  can  be  written  as: 


1 


-DC 

dz 

2 


*  L, 


(23) 


where  Cp,  and  L,  denote  the  specific  heat  of  air  at  constant  pressure  and  the  latent  heat  of 
evaporation,  reflectively. 
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Tenns  1  in  (2.4)  and  (2  J)  denote  die  local  time  rates  of  changes.  As  its  name  in^lies,  the  CLUSE 
model  calcttlates  die  stationary  solntitHi.  Hence,  tenns  1  are  zero  at  equilibtium. 

Terms  2  in  (2.4)  and  (2J)  rqneseiU  die  diffusion  by  die  air  turbulent  movements  and  are 
modelled  by  of  exchange  coefficients  Kj.  The  coefficients  are  equal  for  all  droplet 

categmies,  water  vqxir  and  temperature  and  are  given  by  die  eddy  viscosity  o,: 


3o~ 

PjM'  -  -  K-^  Y  ■  n,v 

^  (2-6) 

az 

K  -  Ot  /  0.74 

where  the  factor  0.74  represents  the  turbulence  Schmidt  (water  vapor,  droplets)  or  Prandtl 
(tenqierature)  numbers  that  are  all  equal  (see  also  [Melville  and  Bray,  1979]). 

Term  3  in  (2.4)  denotes  the  counterdiffiision,  uhich  appears  because  die  droplets  do  not  exactfy 
follow  die  turbulent  motions  of  die  air  doe  to  dieir  inertia.  Tenn  3  reduces  the  turbulent  diffiurion 
as  a  function  of  droplet  radius  and  turbulence  intensity.  It  is  modelled  by  an  exchange  coefficient 

K„': 


Pn^L 


Kn  =  -Kn 


1  - 


1  C,  (V|.n  /  1.56U?) 


(2.7) 


where  Vf  and  u*  dmote  final  fall  velocity  and  friction  velocity,  reflectively.  The  constant  Cj  is 
one  of  the  two  model  tune  parameters  to  be  adjusted  by  comparison  with  experimental  data. 


Tom  4  in  (2.4)  governs  die  macroscopic  fluxes  due  to  ejection  of  dxipl^  from  die  surfiu«  and 
dieir  gravitational  fall  Rouault  etaL  [1991]  proposed  to  model  diis  term  Ity: 
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a<pK> 


^iTfly.n 


[  Pb<*>  -  pS(2)] 


^ere  Tjiy  denotes  die  flight  time  of  die  diqilets  and  die  cracentratioD  in  non-tmbulent  and 
non-evqiorative  conditions.  The  constant  C]  is  die  odier  model  tune  parameter  to  be  adjusted  by 
conqiatison  with  eiqierimental  data. 

Finally,  term  5  in  (2.4)  represents  the  source  terms  due  to  droplet  evaporation.  Widiin  the 
hypothesis  of  cmservation  of  die  total  number  of  droplets,  Mestayer  and  Lefauconnier  [1S>881 
formulated  a  model  diat  yields  an  expression  for  Sb*. 


(2.9) 


In  the  CLUSE  model  for  fresh  water  the  rate  of  change  of  droplet  mass  r  is  described  by  Beard 
and  Plnqipacber's  [1971]  model  of  an  evaporating  droplet  falling  at  terminal  velocity,  which  also 
gives  good  results  for  evqxnating  drqilets  in  a  turbulent  flow  [Edson,  1989]. 


When  we  substitote  (2.6)-(2.8)  in  (2.4)  and  (2.5),  the  set  of  CLUSE  equations  is  closed.  The 
CLUSE  equations  all  have  the  same  general  form: 

.  W  ’  (^■^)  *  S*  -  0  (2 ,0) 

*  Sp  + 

^i^iere  S«  includes  every  term  but  turbokiit  difiiisioo  and  inertia  coonterdifhisioa.  The  stdothm, 
Le.,  die  time  limit  of  dw  non-statkmaty  equation,  is  found  by  an  iterative  process.  The  equations 
are  discretised  by  the  finite  volume  mediod  [Patankar  and  Spalding,  1970]  and  subsequenffy 
solved  widi  a  tridiagonal  algotidun  [SchiesteL  1991]. 


Initial  and  boundary  conditioos  ate  needed  to  solve  die  eqnatioiis.  Fbr  water  vapor  and 
temperature  these  have  been  defined  as; 
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Initial: 

Uppex  boundazy: 
Lower  boundazy: 


Logarithmic  profiles 

“  T(2^)  P^(2^) 

T(W  - 

^•urf 


Pv<Z**) 
Pv«  ^ 


whoe  pv,  denotes  the  satnrated  water  viq>ar  ctmcentratioo.  Fbr  the  droplets  diese  ccmditions  are: 

[  Initial:  pI^z)  >  0 


Upper  boundazy:  p^Cz^J 


Lower  boundazy:  Pa<z,in) 


0  and 


where  and  V^,  denote  the  concentration  of  ejected  droplets  and  die  ejection  velocity, 
respectively.  The  boondaiy  ctmditions  result  firom  die  constraint  diat  die  total  number  of  droplets 
throughout  die  domain  is  conserved.  At  die  lower  boundary,  die  concentration  of  droplets  is  the 
sum  of  the  concentrations  of  droplets  ejected  from  the  sea  and  falling  back  into  die  sea. 
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3  MOIffiL  FOR  THE  AIR  FLOW  OVER  THE  WAVES 

3.1  Vigiulizatk»  of  tbe  calculation 

In  this  duqpter  we  disciiss  our  2-diinensioaal  model  of  the  air  flow  over  die  waves  md  <me 
hmizoatal  coordmale).  The  model  will  be  used  to  cslcnlate  die  hajectorin  of  sea  tpny  thoj^ts 
and  sttbsequendy.  die  vertical  ctmcentiation  profiles  of  the  droplets. 

It  is  in^xH^tant  to  ke^  in  mind  that  die  SeaChise  model  is  an  1-dimensional  Euleiian  model 
(height  only)  and  diat  it  does  not  deal  with  the  coordinates  x  and  y.  This  can  be  visualized  by  tai 
observer  who  remains  fixed  at  a  certain  position  x  and  watches  the  waves  go  by.  The  waves 
induce  changes  in  the  values  of  U(z),  T(z),  p,(z),  etc.  Since  the  model  calculates  a  stationary 
solution,  die  obsover  cannot  specify  diese  variables  as  a  function  of  time.  Therefore,  diey  have  to 
be  averaged  over  the  waves,  Le.,  to  be  averaged  over  a  time  t »  or^  where  oi  is  die  wave  period. 

In  our  model,  variables  like  U(z)  vary  periodically  in  time  with  co  and  it  suffices  to  average  over  t 
-  ar^.  This  is  most  conveniaitly  done  in  a  coordinate  system  that  moves  along  widi  die  wave 
(widi  the  phase  velocity  C  J.  bn  such  a  system,  die  position  of  dw  wave  surface  does  not  move,  or, 
tl(x,t)  =  Tl(x).  The  integratkm  (avenging)  in  time  mentioned  above  is  dien  replaced  by  an 
integration  in  space  (over  die  wave  surface). 

3.2  Basic  parameters 

Our  aim  was  to  develop  a  model  that  only  requires  hqiut  parameters  diat  are  readily  available, 
sudi  as  wind  qieed  and  temperature.  As  a  result,  the  key  input  parameter  for  die  model  is  Uiq, 
Le.,  die  wind  qieed  at  a  referaioe  height  of  10  meters.  Li  our  model,  Uiq  directly  determines  die 
anqilitude  and  frequency  of  die  wave,  as  well  as  the  friction  velocity.  The  latter  parameter  plays 
an  mqioitant  role  in  the  calculation  of  die  vertical  profiles  of  U,  W,  p,  and  T,  as  we  will  see  later. 

The  frktkm  vdodfy  n«  is  calculated  fimn  Uio  and  die  drag  coeffidoit  Qm: 

U?  -  Cj,  Uio  ^3  jj 
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Hm  vahw  of  the  dng  coeffickat  has  been  measmed  many  auOan  (see  [Geenaert,  1990]  for 
an  overview)  and  dqiends  on  Uiq.  wave  stke  and  region.  We  nse  vahies  for  open  ocean  (Notdi 
Atlantic)  and  coastal  regions  (Meelpost  Nbcrdwijk)  as  given  by  [Smidi,  1980]  and  [Davidsm  et 
aL,  1993],  re^ectivety: 

{10^  *  0.61  0.063  open  oce2m 

10^  Cq,,  ■  0.43  •«-  0.100  Ujo  Meetpost  Nooxdwijk 

For  tbe  first  version  of  die  SeaQuse  model,  we  have  chosen  a  single  wave  model  We  assume 
diat  die  waves  can  be  described  by  a  first  order  Stokes  wave: 

H(x,t)  =  -  -l-H,  cos(kx-ut)  +  -l-Hok  cos(2kx-2ot) 


wliere  the  wave  aiqilitode  Hq  and  fiequency  oi  are  given  by  [Edam,  1990]: 


u 


g 

riu„ 


Ho 


*  ^  Pio 

log 


(3.4) 


The  factor  Fj  in  (3.4)  is  die  wave  age  parameter,  which  in  our  model  is  fixed  at  F]  ■  0.7,  Le,  die 
vahie  for  fully  developed  waves.  Obviously,  die  assumption  of  fully  developed  waves  will  often 
be  violated  in  coastal  tegkms  (see,  e.g.,  [Maat  et  al,  1991]). 

3.3  The  eddy  difiusivity 

The  eddy  diBaavity  is  Inportant  parameter,  since  it  governs  die  difTbsion  of  heat,  moisture 
and  droplets  by  die  air  turbulmt  movements  in  die  CLUSE  model  (see  equathm  Q.6)).  In  additi¬ 
on,  the  eddy  diffosivity  qipears  indirtcL’"  in  ^7)  for  die  counteidifiusion  due  to  droplet  inertia. 

To  calculate  die  eddy  diffosivity,  we  divided  die  vertical  axis  in  two  regions: 

(a)  The  wave  region,  vhere  die  waves  influence  die  air  flow. 

(b)  The  atmoqiheric  region,  vdiere  die  waves  no  longer  'nflnecce  me  air  flow. 
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In  die  atmcMpbcfic  r^kxi.  we  apply  die  condition  of  a  folly  devekqied  boundaty  layer  (classical 
atmospheric  sorfiax  %er).  b  diat  case,  die  eddy  difiiissviQr  0(  u  gtven 

0((z)  -  K  u.  z  z  >  Ho  (3J) 

hi  die  wave  r^ioo,  we  aasonie  diat  momentum  is  «ity  transftned  to  die  water  and  not  to  the 
waves.  We  also  aasmne  diat  the  stractme  at  die  mean  flow  is  localty  diat  of  a  turbulent  boundaty 
layer  over  a  flat  surface.  In  diat  case,  local  profiles  afitjt)  may  be  calculated  by: 


Ot  (x,  z)  ■  K  u,  (z  -  H  (x) )  n  (x)  <  z  <  Ho 

Oo(x,z)  -  0  --l-Ho  <  z  <  ii(x) 


(3.6) 


vriiere  Ti(x)  coneqxmds  to  die  bei^t  of  die  wave  surface  at  position  x.  Note  diat  T|(x)  is  not  a 
function  of  tinie,  since  we  perform  die  cakulatkn  in  the  coordinate  system  diat  moves  along  with 
the  waves.  The  average  profile  o,(z)  in  die  wave  region  is  calculated  by: 


Ot(z) 


0((x,z)  dx 


-|Ho  <  z  <  Ho 


(3.7) 


vdiere  it  suffices  to  integrate  over  half  a  wave  period,  Le.,  from  trough  to  crest,  due  to 
symmetry.!^ 

h  has  been  rqiotted  diat  die  wave  regkxi.  uiieBe  die  waves  influence  die  air  flow,  extends  q>  to  5- 
10  times  H,  (Ota  wave  crest  is  at  z  ■  IdHo).  However,  when  constmctiiig  die  profile  of  a,  by  (3.5) 
nd  (3.6),  we  found  diat  diere  was  no  difference  between  die  average  of  die  local  profiles  and  the 
'atmoqpheric' Ot  above  z  ■  Hq.  Rgure  3.1  diows  the  resulting  profile  of  o,(z).  The  net  effect  die 
wave  surfice  int^iation  is  diat  die  slope  of  a/z)  is  decreased  and  consequent^,  o^z)  *  0  is 
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Uat  k  aqaidiMadjr  weead  «*w  te  wave  Mrfhoe,  wbidl  k  aot  We  ana  at  a  fM  efaldiittailjr  ^aeed  la  X.  ta  wpciiilk  B  we  Aaw  te  We 
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readied  below  z  *  0.  The  modier  of  local  profiles  dial  are  avenged  is  not  very  inqxxtant  A 
mmifinini  of  20  is  required;  averaging  500  profiles  yields  essentially  die  sanaeiesiilts.Pl 


ngui«3.1:  Veiticalprofik  of  eddy  diCbuivity  in  die  loweriegkn  of  the  cakalatioa  domain 


3.4  The  profile  of  horizontal  wind  qieed 

With  die  hypodiesis  that  the  oonstant  flux  assunqition  is  valul  for  nwwnwifiim  in  die  atmospheric 
region,  die  wind  profile  be  calculated  from  q,(z)  by; 


z  >  Ho 


(3.8) 


adiere  ^  denotes  die  atmoqdieric  stabilfry  funcdoo  (see,  e.g.,  [Edson  et  aL,  1991])-  In  neutral 
stratification  4b  **  U  ^  noii*neatral  oonditkms,  die  Monin-Obhnkov  lengdi  L  is  needed  to 


AfpMaia  C  dhcaMM  anote  MM  to  cdodM*  S«e<  aMB^pdoa  or  a  lUtjr  SmSoped  «S  doaiiual  wave  mor.  Hmww.  wtt 
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calculate  4b-  Becanae  zq,,  T*  and  (V  SR  not  3^  availaUe,  it  is  not  straigjitfbrwaid  to  obtain  L. 
Therefore,  SeaChiae  is  developed  fiat  for  neutral  conditions. 


Equation  (3.8)  shows  that  a  reference  value  for  U  is  needed  to  constiuct  die  profile.  For  this,  we 
take  die  value  of  the  hqmt  variable  Uio  at  z  b  10  m.  Rom  diere,  the  wind  profile  is  constructed  to 
die  top  of  die  boundary  layer  and  down  to  z  *  Hq. 


Tlw  wind  profile  dins  obtained  can  be  used  to  infer  anodier  inqxntant  parameter  for  the  nindellin£ 
of  tufbulence:  die  roughness  length  Zq..  To  diis  end,  linear  regression  is  ^iplied  to  fit  die  wind 
profile  to: 


o(z).^ln(.^)  z>H. 


(3.9) 


which  yields  a  value  for  Zq,.  We  may  conqiare  die  value  of  Zq,  widi  values  for  die  roughness 
kngdi  Zq  obtained  by  Chamock's  relation  or  die  Kitaigoroddd  model  (see  qipeadix  D).  Up  to 
now,  we  have  observed  dut  Zo,  is  nearty  equal  to 

The  wind  profile  in  die  wave  region  is  calculated  in  die  moving  coordinate  system,  vdiere  T)(x,t)  = 
T|(x).  in  diat  case,  local  logarithmic  profiles  U(x,z)  are  given  by: 


U(x,z)  -  Uoi,(x) 


In  ~  ^tx)  Zq^Cx) 

K  \  Zot(x) 


n  (x)  <  z  <  Ho 


(3.10) 


OolW  denotes  die  local  value  of  U(x,z)  at  die  wave  sutfine  and  Zq|J(x)  represents  the  local' 
roughness  lengdL  It  is  easily  verified  diat  u^f*)  ■  U(x,t)(x))  when  z  -  ti(x).  The  wind  qieed  at 
die  surface  is  composed  of  a  local  component  due  to  die  water  orbital  movement  under  die  wave 
and  a  constant  oonqionent  due  to  suifiice  sl^: 


U(x,i)(x))  -  U1)(X)  ♦XU, 


(3.11) 
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(D  it  die  fkeqaeacy  of  die  Stokes  wsve.  Sabseqneiidy,  we  apf\y  the  oonditioo  that  all  local 
profiles  U(x^)  most  converge  to  the  vafaie  of  Uiq  at  z  ■  10  m,  Le.,  outside  the  wave  region. 
InqxMing  diis  condition,  ^(x)  can  be  calculated  by: 


^10  “  Uoi.(x) 


U.  flO  -  n<x)  ♦  2oi.(x) 

K  \  *01.  (x) 


(3.12) 


Finally,  die  average  profile  U(z)  in  die  wave  region  is  calculated  by: 


u<2)  .//"ocx.z)  dx  -Ah.  <2<h. 

vdiere  dx  tuns  over  half  a  wave  period,  i.e.,  from  trough  to  crest  The  number  of  local  profiles 
should  be  fairly  large  (»  5(X))  to  obtain  a  smooth  wind  profile  in  the  wave  region.  Figures  3.2a-c 
show  T)(x),  Ool(x)  and  Zol(x)  for  Ujo  «  10  m/s.  The  largest  values  of  ox  and  zql  found  at  die 
crest  Note  that  die  mean  value  of  u^^  slighdy  larger  than  Ku*,  due  to  die  asymmetry  of  the 
Stokes  W8ve.(^l 

Rnally,  figure  3.3  shows  the  conqilete  wind  profile  U(z)  for  Uiq  »  10  m/s.  A  nearly  perfect  ™trli 
is  obsoved  at  the  intersection  of  the  atmospheric  and  wave  regions  at  z  »  Hq. 
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The  teo^Mntnre  prafik 

llie  initial  profile  of  temperature  is  cakalated  in  almost  die  same  manner  as  die  wind  profile. 
Thus,  in  die  atmospheric  region,  where  the  constant  fhn  assomptioo  is  valid  in  die  absence  of 
evaporating  droplets,  a^t)  is  oaed  to  calculate  T(z): 


dr  _  U.  T.  Pie 

*3? - F, — 


(Z  >  Ho) 


(3.14) 


sriiere  die  ftandd  number  IV,  «  0.74  [Rouault  etal.,  1991].  The  model  input  parameto-  Tjq 
(tenqierature  at  10  m  AMSL)  serves  as  the  reference  value  needed  to  calculate  the  profile.  The 
scaling  factor  T*  is  inftned  fmru 


T,  ■  ^  (Uio  ■“  Z'i*)  (Tjo  ~  ^auzf^  (3.15) 

where  we  use  die  average  value  »  0.001. 

The  tenqierature  profile  in  the  wave  regitm  is  again  calculated  in  the  moving  cocndinate  system, 
from  local  logarithmic  profiles  T(x,z): 


T(x,z)  - 


Pit  T,  fz  -  ii(x)  »  Zo,.(x) 

K  \  Zoi.(x) 


n  (x)  s  z  <  Hp 


(3.16) 


uhere  we  assume  that  die  surface  tenqierature  does  not  vary  over  die  wave.  The  local  tenqierature 
roughness  lengdi  Z()l(x)  does  dqpend  on  x  and  is  again  found  by  inclosing  the  condition  that  all 
local  profiles  converge  to  T|o  at  z«  10m: 


Dr 


*10 


(3.17) 


The  average  profile  T(z)  in  die  wave  tegkn  is  calculated  by  integrating  the  local  profiles  over  die 
wave  surfine.  Just  as  m  die  case  of  die  wind  profile,  the  number  of  local  profiles  riiould  be  circa 
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500.  Hgure  3.4  shows  die  temperature  pnifik  for  Tio  «  20  *C  and  «  10  *C  (and  Uio  =  10 
m/s).  Note  that  die  atmoqiheric  profile  and  the  wave  legkni  profile  match  perfectly  at  die 
intersectum  (z  »  Hq). 


10.0  tzs  ISO  ira  20a 


ingttre3.4:  FtofileofteoqieraturBasafunctiaat^beigfat 


3.6  The  humidity  profile 

The  humidity  profile  of  humidity  is  calculated  in  exaedy  the  same  manner  as  the  ten^ierature 
profile.  Thus,  for  the  atmoqihaic  region; 


3pv  .  u.  Pv,  SCt 

~Sz 


(Z  >  Hp) 


(3.18) 


where  die  Schmidt  number  Sc, »  0.74  [Rouault  etaL,  1991].  The  scaling  factor  p^*  is  calculated 
from: 
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Pv. 


(U„  -  KU.)  (p^  „  -  p^  „) 


(3.19) 


we  use  die  average  value  Cb  »  0.0012.  The  values  fin:  pv.io  aixl  Pv.o  ^  obtained  fran 
conversion  of  the  iiqiut  parameter  RH|o  and  the  relative  humidity  at  die  sea  surface  (RH  ••  97%, 
determined  by  the  salinity),  reqpectivefy. 

In  the  wave  region,  local  logarithmic  profiles  Pv(x,z)  are  calculated  with  q  s  RH,aif  ^ 
roughness  length  Z(b,(x)  inferred  from  py  lo  by  expressions  analogous  to  (3.16)  and  (3.17).  The 
average  profile  pv(z)  is  calculated  by  an  integration  over  circa  SOO  local  profiles  (efr  (3.13)). 

An  exanqile  of  the  resulting  overall  profile  of  humidity  is  shown  in  figure  3  J,  for  RHio  =  75% 
(and  Uio  »  10  m/s,  Tjq  =  20  *C,  Tq  «  10  'C).  Fbr  cmvenience,  the  profile  is  given  as  a  relative 
humidity  profile. 

6acs707s«oas«ogsioo 


Figare33:  Profile  of  idadvefannudity  as  a  function  of  bei^ 
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Figure  3.6:  Fanuneten  for  the  wind  field  W{x^) 

a)  Height  of  wave  surface  versus  x>giid  index 

b)  Piston  velocity  versus  x-grid  indn 

c)  Example  of  local  profile  of  vertical  wind  component 

d)  Vertical  wind  con^xnent  distriboticn  over  foe  wave  at  10  m  (see  text) 


Figure  3.6  shows  results  for  a  cakulation  of  W(x^).  In  fliis  specific  case,  Uw  =  10  ni/s  and  1000 
local  U-profiles  were  used  (x-grid  density).  The  wave  travels  in  the  +x  direction  and  foe  z-grid 
has  Nz  grid  points  between  z  =  z  *  10  mtr.  Figures  3.6a  and  3.6b  foow  foe  height  of  foe 

wave  surface  ti(x)  and  foe  pisttxi  velocity  as  a  function  of  x-grid  index,  reflectively.  Hgure  3.6c 
shows  a  local  wind  profile  W(x,z)  at  foe  windward  side  of  foe  wave  (i.e.,  wave  surface  moving 
down),  obtained  with  NzslOCX). 
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The  success  of  die  cakulation  of  the  locsl  jsi^les  W(x^)  may  be  verified  in  die  region  well 
above  die  waves,  where  die  waves  no  longer  have  any  influence  cm  W(x,z).  Tims,  in  that  r^ion 
W(x,z)  should  be  independent  of  x.  Rgure  3.6d  shows  W(x,z)  at  a  height  of  10  mtr  (■>  6  J  Hq  for 
U|(^10  nt/s)  as  a  function  of  x.  It  was  constructed  by  taking  the  value  W(x,z*10)  of  each  of  die 
local  profiles  W(x,z).  The  curve  by  was  obtained  firom  local  profiles  W(x,z)  diat  had 

been  calculated  widi  1000  vertical  grid  points  between  the  surface  and  10  m  AMSL,  the  curve 
marked  by  (x)  from  local  profiles  W(x,z)  diat  had  10000  vertical  grid  points.  It  is  easity  seen  from 
figure  3.6d  diat  Nz  should  be  in  die  order  ctf  10*  to  obtain  a  wind  field  W(x,z)  diat  is  independent 
of  X  above  the  wave  region. 


Since  bodi  U(x,z)  and  W(x,z)  are  now  determined,  a  vector  plot  of  die  wind  field  over  die  wave 
can  be  constructed.  Such  a  plot  is  diown  in  figure  3.7,  where  die  length  and  direction  of  the 
arrows  rqiresent  wind  speed  and  wirul  direction,  respectively.  Our  model  for  die  air  flow  over  die 
waves  results  in  a  peifocdy  symmetric  wind  field  with  a  rotor  between  die  wave  crests.  The 
influence  of  the  waves  decreases  widi  height  and  becomes  negligible  above  2-3Ho. 


Rgure  3.7:  Air  flow  over  die  wave 


TNOraiMrt 


nige 

30 


4  DROPLET  MOVEMENT  IN  THE  AIR  FLOW 

4.1  Ejectkn  hei^  and  ejecticm  vekxd^ 

Li  order  to  calculate  die  trajectories  of  die  sea  spny  droplets  diat  result  from  bursting  bubbles,  we 
need  to  know  their  ejection  velocities  V^(D).  To  diis  end,  we  first  infer  die  ejection  heij^ 
H^j(D)  of  die  driblets  from  Blanchard's  ejection  data  [Blanchard,  1963,  1989]  for  jet  droplets 
consisting  of  distilled  water  of  22-26  ‘C.  Blanchard's  figures  show  negligibte  difference  between 
distilled  water  and  seawater.  Snbsequendy,  ejection  velocities  V|^(D)  are  r-almlatMt  from  Hg<D) 
by  an  iterative  procedure  described  previously  [Wu,  1979;  Rouault  etal.,  1991].  The  method  uses 
the  equation  of  vertical  motion  in  non-turbulent  and  ncm-eviqiorative  conditions  over  a  horizontal 
still  surface  (see  below,  (4.2))  and  an  initial  guess  for  V|^(D)  to  infer  the  marimum  height  H(D) 
that  the  droplet  reaches.  The  value  for  V^(D)  is  dien  gradually  adjusted  untill  H(D) »  Hg<D). 

4.2  Drt^let  trajectories 

The  presence  of  waves  complicates  die  calculation  of  die  drcplet  trajectories,  because  the 

horizontal  displacement  of  die  dixplet  and  the  movement  of  the  air  itself  most  also  be  taken  into 
account  The  trajectories  are  most  conveniently  calculated  in  the  coordinate  system  that  moves 
along  with  die  waves  (phase  velocity  C,).  hi  diis  tystem,  die  horizontal  and  votical  droplet 
velocities  o^  and  w^  are  given  Ity: 

{Uj(D,x,z,t)  «  u,(D,t)  +U(x,z)  -  C, 

Wi(D,x,Z,t)  -WjfD.t)  +W(x,z) 


where  U(x,z)  and  W(x,z)  are  die  local  horizontal  and  vertical  wind  qieed  (see  §3.4  and  §3.7)  and 
u^,t)  and  Wf(D,t)  denote  the  horizontal  and  vertical  velocity  of  the  drtplet  relative  til  the  air. 
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ia  (4.2),  Q>  denotes  die  drag  coefiBdent;  for  a  qdierical  drtqdet  is  a  fdoctioa  of  die  Reynolds 
mimber  Re: 


Re  =  *  w’  ^ 


(4.3) 


ti^ieie  V  is  die  kmemadc  viscosity  of  air.  The  drag  coefficient  is  \rf  [Raudvild,  1979]: 


Co  -  24  /  Re 

Co  «  24  (1  0.19Re}  /  Re 

Co  -  24  (1  +  O.lSRe®  *”)  /  Re 


Re  <  0.5 
0.5  <  Re  <  2 
2  ^  Re 


(4.4) 


The  initial  velocities  u,(D,tsO)  and  w^,tsB0)  needed  to  calculate  u^.t)  and  w^,t)  by  (4.2)  are 
given  by: 


I  u,(D,t-0)  -U(x,j,n(x.3))  «Uoo(x.j) 
I  w,(D,t=0)  -V^(D) 


Note  that  the  initial  velocity  1^(0, t=0)  depends  cm  the  position  over  die  wave  surface  where  the 
droplet  is  ejected  (and  is  equal  to  the  horizontal  conqionent  of  die  air  flow  at  the  surface). 

Equations  (4.1)-(4.S)  allow  die  calculation  cS  die  trajectory  of  a  droplet  ejected  at  x  s  The 
calculation  is  aborted  when  die  droplet  hits  die  water  surface,  i.e.,  iKlien: 

2d  ^  n  (Xd)  (4.6) 

The  calculation  is  most  effectively  done  in  a  grid  diat  is  equidistantly  placed  along  die  droplet 
trajectory  .1^  Such  a  grid  widi  stqi-size  Sp(D)  is  built  recalculating  At  in  each  cycle: 

At(D,x,z,t)  ■  —  ftp(D)  .  — ^ 

^  Uj(D,X,Z,  t)*  +  W4(D,X,Z, t)’ 


The  B^eeUfy  on  ako  be  ctleabMcd  ia  the  mace  oonBOB  mwner.  Smi  ii,  widi  a  coosmil  at  (tniiadly,  at  >  so  |U  a  needed),  bM  in  dw  caw 

S-10  JmetaamanyM^me  imBifCd  lO  Olittfa  the  SMBC  iccwiKy. 
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The  rimtt  ttep  At  u  «m»n  when  die  velocity  of  die  dro]^  with  reject  to  die  wave  nnface  is 
UghJ^  s^Mreas  it  becomes  larger  s^ien  die  drqiilet  slows  down.  The  stq>  size  8p(D)  depends  on 
the  size  of  die  dn^let,  and  varies  between  0.02-0  J  cm  for  20  ^  D  ^  SOO  tun. 

An  exanqile  of  a  typical  trajectory  is  shown  in  figure  4.1a.  In  diis  case,  Uio  ^  10  m/s,  die  wave  is 
travelling  in  the  -t-x  direction  and  a  dn^Iet  widi  a  /Kamnrr  of  10'  'm  was  ejected  firom  die  wave 
trough.  It  is  anjdiasized  that  this  trajechny  is  calculated  in  a  nou-turbulent  and  non-ev^torative 
atmoqihere. 

Different  trajectories  are  obtained  when  droplets  (D  ^  1 10  pm)  are  ejected  fiom  the  tee  side  of  die 
wave.  They  do  not  siiiqily  fall  back,  but  are  lifted  iqiwards  and  enter  in  die  wave  rotor  (figure 
3.7).  In  some  cases,  diey  return  to  the  surface  (see  figure  4.1b),  but  in  odier  cases  th^  start  to 
move  in  a  circle  (see  figure  4.1c). 


« 


lUi  DMiii  dm  the  ralMhc  vaiociqr  or  the  draplei  a  hi|h.  (be  droplet  k  M  a  potWoD  (x,z)  where  the  wind  fidd  (U,W)  to  ittMW. 
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ngiiie4.1:  Represeattiive  droplet  tmjecto^ 

«)  tt:^ectaty  iMt  influenoed  by  dK  wive  rotor 

b)  Thijectaiy  entering  the  wive  rotor 

c)  Giculir  tnjectoiy 


Gearly,  die  aitenmn  that  the  trajectoiy  calculation  is  finished  when  die  droplet  hits  die  siiifiK« 
again,  will  not  work  for  droplets  moving  is  a  closed  tnyectory.  In  onr  opinion,  die  circular 
trajectories  are  an  artefact  resulting  from  the  perfect  symmetry  of  our  model  of  die  air  flow  over 
die  waves.  In  reality,  the  wave  surface  is  not  perfecdy  smoodi.  Small  variations  in  Ti(x)  induce 
aqmmnetry  in  the  wave  rotor  with  die  consequence  diat  the  drtqilets  do  not  move  indefinitively  in 
a  perfect  circle.  Therefore,  we  decided  to  stop  the  calculation  of  circular  trajectories  after  a  full 
circle,!^  widi  die  simultaneous  removal  of  die  droplets  from  die  atmoqihete. 

4.3  Flight  time 

The  suspeosko  or  flight  time  of  die  droplet  is  die  sum  of  all  dme  intervals  At  of  die 

trajectory  cakulatwHi.  A  {dot  of  as  a  frmction  of  die  position  x  where  the  droplet  was 

ejected,  is  shown  in  figure  4.2a.  In  diis  case,  Uk^IO  m/s,  I>bS00  pm,  x^l  and  xw51  coneqiood  to 
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wave  trou^  «iH  cieat  leipecUve^,  and  die  wave  is  travelling  in  the  ■I'X  direction.  The  fli^t  time 
is  larger  at  the  lee  side  of  die  wave  (x  >  SI,  snrfKe  moves  up)  due  to  die  airflow  thtt  moves  the 
air  upwards.  The  opposite  effect  is  observed  at  die  windward  side  of  die  wave,  where  the  air 
moves  down  and  dmlens  die  flig^  dme. 


Rgare4.2:  Fl^  times  T*.  over  the  wave  (lee  text  for  emanation  of  discoiitinuities) 

a)  ^  for  dropm  of  DbSOO  |un 

b)  T^  for  dro{4ets  of  I>il20  |un 

c)  for  drofdets  of  !>■  20  |im 

Figure  4.2b  shows  Tq^xJ))  as  a  function  of  x  for  a  droplet  of  D^IZO  pm.  In  diis  case,  the 
prolongatkio  of  the  flight  time  at  die  lee  side  of  die  wave  is  much  larger  dian  the  shortening  m  die 
weadier  side.  This  reflects  die  facts  diat  die  droplets  ejected  at  tibe  weadier  side  eater  die  wave 
rotor  and  conqilete  a  semi-drcnlar  trajectory  (see  figure  4.1b). 

As  die  diameter  of  die  droidet  decreases,  die  prolongation  of  die  flight  time  becomes  larger  and 
larger.  Rgure  4.2c  diows  Tq^xJ))  for  droplets  d  I>*20  jim.  The  general  shape  of  the  cove  in 
fig.4.2c  can  be  eiqilaiiied  in  terms  of  die  wave  rotor,  bm  diete  are  also  two  artificial 
discontinuities  in  fig.4Jc  («  je-SO  and  xw68)  due  to  the  discretization  of  the  trajectory 
calculation. 
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The  x-dq>endency  of  J))  is  nmoved  in  die  integintiati  over  die  wave  amhoe: 

Tfiy(D)  dx  (4*) 

A  plot  of  T||y(D)  vemis  D  is  shown  in  figme  4.3  (carve  madced  by  (*)),  vdiete  we  also  present 
T||^)  for  the  case  of  a  flix  smfoce  (curve  marked  by  (x)).  h  is  easily  seen  diat  the  presence  of 
waves  has  litde  infhience  on  the  flight  dme  of  the  large  and  heavy  droplett  diat  foil  down  quickly 
due  to  gravitation,  but  bas  a  large  effect  on  die  times  of  die  small  and  li^t  dtr^lets  diat 
may  enter  the  wave  rotor  and  start  to  move  in  a  circle. 


0  MO  aoo  M»  400  aoo 


Figate4.3:  Average  flight  times  for  dnqileis  over  the  waves  (+)  and  over  a  flat  ■arfooe(x) 

The  number  of  x-grid  points  over  die  wave  (Nx)  needed  to  tditain  die  average  flight  T0y(D) 
nmges  fiom  only  20  in  die  case  of  large  droidetB  to  50  for  droidets  dud  taKy  eater  die  wave  rotor 
(D  <  100  pm),  hi  foct,  a  fordier  increase  t^Nx  does  not  yield  more  reliable  average  flight  times. 
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4.4  Rnal  fidl  vdodties 

The  pkMs  of  V|(xJ>)  vemu  x  are  leas  cooiidicated,  tmce  aH  exoqit  die  largest  dni(dets  readi 
final  faU  vdodty  at  a  ladier  eariy  stage  of  dieir  tnyectosy.  Therefore,  V((xJ>)  is  only  dqwodeat 
on  X  for  die  largest  droplets  (see  figure  4.4  for  a  plot  of  V|(xJ>)  versos  x  for  droidets  of  D«SOO 
pm).  Not  smpriaingty.  the  larger  values  for  V|(x4>)  are  found  at  die  lee  side  at  the  wave  where 
die  flight  times  are  koger. 

Hgure  4  J  shows  die  average  final  fdl  velocities  V|(D)  versus  D  after  integration  over  die  wave 
surface  (+ )  (Nx  as  for  die  fli^  times,  see  above).  For  compariscm,  die  ftU  velocities  in  die  case 
of  a  flat  surface  are  also  {dotted  (x).  The  two  carves  are  nearly  identicaL 
0  ieaoao40Soso7DMioiw 


Figiire4.4:  AialCdlvelocities  Vf  over  the  wave  for  droidets  of  D-500  pm 
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FSgure4^:  Avenge  final  &D  velocities  for  droplets  over  the  wivea(-f)  and  over  a  flat  amfiKcCx) 

4.5  I^ofiles  of  droplet  amcentratiao 

The  last  and  most  conylicated  of  the  dn^det  parameters  is  the  vertical  profile  of  driblet 
concentration.  This  profile  is  calculated  with  our  model  of  die  airflow  over  die  waves,  and 
therefore  iqMesents  noo-tutbuleat  and  ntm-ev^xnative  conditions.  It  is  used  in  SeaOuse  as  a 
relaxation  term  in  die  modelling  of  the  macrosct^  flux  of  die  droplets  due  to  ejectkm  and 
gravitatkaial  fidl  (see  (2.8)  in  chqiter  2). 

Rrst,  let  ns  consider  a  sin^  droplet,  ejected  at  a  positkn  X  at  the  wave  surface.  Hie  probabili^  ^ 
of  finding  the  droplet  in  die  beig^  interval  [z,z4Az]  is  related  to  die  time  T(z,z-»At)  diat  die 
droplet  qiends  in  diat  height  interval: 


T(z,Z4-Az 

Tny<D) 


(4.9) 


TNOraport 


nige 

38 


I 


The  x-ckpendeacy  is  removed  by  taking  die  anm  of  the  coattbatiaiis  of  individiisl  droplets  ejected 
at  various  poiitiofis  X  over  die  wave  tmfKe: 


Cd(2,z+Az)  -  Cd.*(*»**Az) 

Which  is  renonnalized  by: 


(4.10) 


Cd<2) 


lim  Co(z,z*Az) 
Az^O  # 

J  Cd(2»*+Az)  dz 


(4.11) 


The  dropl^  cooceotration  profile  po^z)  r .  non-evqxirative  and  non-turbukiit  conditions  is  found 
by  mult^lkation  widi  the  density  Pej(.D)  of  dn^lets  ejected  at  the  wave  surface: 


P0(2)  ■P.j(D)  Cd(*)  (4.12) 

The  density  p«j(D)  win  be  discussed  in  cluqiter  6.  At  diis  stage,  we  restrict  the  cakulation  to  the 
probabilities  (^(z). 

The  interval  Az  in  (4.9)  is  fixed  to  die  maximom  stqvsize  Sp(D)  in  die  trajectory  calculation  (see 
§4.2).  This  dioke  ensures  that  die  dropl^  never  moves  more  dian  one  height  interval  Az  during 
one  st^  At  of  die  trajectory  calculation.  This  facilitates  die  computation  of  Cd(^)>  ^  ^  ^ 

>  qiecific  reason  vdiyAz  cannot  be  smaller  dian  the  maximum  value  of  £p(D). 

The  number  of  individual  profiles  diat  are  averaged  to  remove  die  x-dqiendency  (see  (4.10))  wa8 
first  set  to  ^  (50  for  D<100  pm).  Figure  4.6a  shows  die  resulting  profile  4d(z)  far  1>>200  pm.  In 
diiscaae,U|o*  10  m/liaiui  die  wave  trough  and  crest  are  located  ittz« -0.73  m  and  z>  40.84  m, 
retpectivety.  It  is  easify  verified  diat  20  individual  profiles  are  not  snfficient  to  obtain  a  smoodi 
average  profile.  In  fitet,  die  individnal  profiles  can  easily  be  recognized  in  figure  4.6a. 

I 

A  radier  large  number  of  individual  ptt^les  is  needed  before  the  contributions  of  die  individual 
i  profiles  are  no  longer  eaaity  recognized  in  die  average  profile.  Even  diea,  the  average  pn^le  is 

i 

L 


p 
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not  veiy  smooth  and  a  Ganssian  smnodiing  rontinel*]  had  to  be  apfdied.  Rgure  4.6b  diows  l^z) 
IV200  and  100  indtvidnal  profiles,  and  after  smoodung. 


Rgure  4.6:  Nonnalized  diQ{det  profiles  after  avenging.  Maxinui  are  probaMy  due  to  die  dioice  ot  the 
wave  profile 

a)  avenge  20  individual  profiles 

b)  avenge  of  100  individual  profiles  and  additional  smoodiing 


The  number  Nx  of  individual  {vofiles  needed  dqpends  on  die  horizontal  wind  speed  Uig  (vdiidi 
determines  the  wave  hei^t)  and  varies  between  SO  (U|o  <  5  m/s)  and  200  (Ujo  >  16  m/s).  In  fiax, 
Nx>b200  is  onty  just  sufiBcieiit  for  dre  largest  wind  qieeds  for  udiidi  die  SeaChise  modd  is 
develqied  (Uiq  •  25  m/s).  However,  widi  Nx«200  and  U,o  >  2S  m/s.  die  CPU-time  needed  to 
cakulate  the  profiles  ^9(z)  for  dropleU  widi  D  <  100  pm  is  very  hi^  (up  to  30-40  min  CFU- 
tnde  for  0*20  pm).  An  average  over  even  more  individoal  profiles  would  result  in  unaccqitable 
computational  demands. 


*  h  ek  witaa.  reeb  SsiMotM  aHl  h  unoo^td  wn  fte  ilwmii  toi  .gp**]- 

affl  -  atMMM]  ~ -f  oOHaR]  „ -I- iriMM 


a[j] 


.Jiiil 


3« 


•1-4 
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5  MOIffiL  FOR  THE  EVAP(»ATK»fC»>  SEA  SPRAY  IHICX>LETS 


3.1  Puent  and  offqxmg  categories 

As  pointed  out  Mestayer  [1990],  salt  water  droplets  must  be  described  by  at  least  two 
parameters:  their  radius  r,  and  a  parameter  describing  dieir  chemical  content.  In  die  SeaChise 
model  dus  is  done  1^  defining  parent  cat^ories,  whkh  contain  the  droplets  widi  ejectkm  radius 
r^i)  ±  Arq<i).  Subsequently,  oCfiqving  categories  widi  radii  rg0(ij)  ±  are  defined,  ^iriiich 

contain  all  drtqdets  whose  current  radius  is  r^ij)  and  whose  ejection  radius  was  r^i).  This  leads 
to  a  model  widi  variables  diat  may  have  iqi  to  3  indexes:  parent  categ(»y,  ofEqiring  category  and 
height  Thus,  V-V(ij,z). 

5.2  Andreas'  model  of  ev^iorating  sea  spray  dro(dets 

Andreas  [1989]  started  with  eq.(13*26)  of  Rnqipacber  and  ICto  [1978]  and  showed  that  for  sea 
spray  droplets  diis  equation  could  be  rewritten  as  a  set  of  duee  equations: 


2  gf  _  Y  ♦  g>. 

R  T(z)  p,,  X  m  -  m. 


(5.2) 


1  ♦  m/niy 
1  ♦  IT 

•  M. 


vhete  p  denotes  die  dnsity  of  die  droplet;  m,  m,  and  m,,  die  total  mass  of  die  dn^det  die  mass  of 
water  and  die  mass  of  salt  in  die  droplet,  respectively;  p^  p,  and  p^  the  densi^  ot  water,  water 
vqmr  and  saturated  water  vapor  density,  leyectivety;  L,  die  latent  heat  of  evaporation;  and 
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M,  die  molecular  weig^  of  water  and  NaCl,  reflectively;  y  die  nionber  of  kms  into  udiich  NaO 
and  [Rnqipaclier  and  Klett,  1978],  [Milkto,  1972]: 

v.(T)  -  v,<,(T)  ♦  S;{T>  yc  (5.4) 


In  a^,(|t/M,) 
V.  |i/M, 


(53) 


Dv(T) 


D^(T) 


r  ^  D,(T) 

2  x  My 

r  +  X  «c  ^ 

R  T 

(5.6) 


ki(T) 


k.(T) 


r  _  k.(T) 

2  It 

[  I  +  A,  X  ttf  p.  Cp.>| 

R  T 

(5.7) 


In  (3.4)-(5.7)  V,  is  die  fipaimt  molal  volume,  ^  tbe  practical  osmotic  constant,  Dy  the  modified 
molecular  difiiisivity  of  water  vapor  in  air  and  k,'  the  modified  thomal  conductivity  of  air.  (Tther 
symbols  are  explained  in  fipeodix  A  and  in  die  text  below. 

5.3  Mestayei's  sinqilified  model  fm  sea  firay  droplet  evfxnation 

Mestayer  [1990]  showed  that  Andreas'  equations  may  be  sinqilified  widi  five  assumptions,  valid 
for  r  >  1  pm,  0  ^  T  ^  30*C  and  RH  ^  75% : 

>1:  The  fiparent  molal  volume  v.  (eq3.4)  is  not  a  function  of  tenfierature.  In  that  case,  die  duid 
equation  of  Andreas  may  be  written  as: 


p(x) 


1  ♦  p(r) 


^(rj_y.  I  „  ,  ^  10-»  pM  sjzl  j 


1  ♦ 


(5.8) 
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where 


£•  .£4^  =  (|i(r)-i  +  l)-i  (5.9) 

P 

S(r)  is  the  salinity  of  the  clnq)let.  Equation  (5.8)  can  be  solved  by  an  iterative  procedure  for  a 
given  value  of  r^ 

>2:  The  practical  osmotic  coefficient  ^  (eqii  J)  is  not  a  function  of  solution  molality  p.  /  M,,  but 
a  constant 

>3:  The  first  term  in  (5.2)  may  be  neglected.  This  introduces  only  1%  error  in  the  teiiqreratuie 
and  humidity  range  considered  here.  Assurrqrtions  2  and  3  allow  us  to  rewrite  the  second 
equation  of  Andreas  as: 

y(r)  =  -  0.6776  |i(r)  (510) 

where  the  values  for  the  constants  have  been  entered  (with  ’^  =  1.1,  the  value  at  equilibrium 
atRH»85%). 

>4:  The  molecular  diffiisivity  of  water  yapot  in  air  D,  (eq3.7)  is  not  a  function  of  temperature. 
In  that  case,  (5.6)  reduces  to: 


p(x) 

S(r) 


m, 

m 


D'(r) 


1 


1 

+  8.10-* 


2.26  10-* 


— -  +4.36  10-* 
r '' 


r-i 


(5.11) 


where,  once  again,  the  values  for  the  constants  have  bear  entered  (widi  T  s  283  K,  and 
D,(283  K) »  2.26  lO-*  m^/k). 
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>5:  Tbe  ri***™!  conductivity  air  k.  (e4*S.8)  is  not  a  functko  of  tempendnre.  bi  diat  case,  (S.8) 

may  be  written  as: 


^  \  J  ^  0.95  z>l|tm 

1-  1.00  r  >  5  im  (512) 


Fbr  tbe  range  of  radii  considemd  here  (r>  1  pm),  we  may  approximate; 
k.  «  k.  -  0.0249  (T  «  283  K) 


(5.13) 


Assunqitioas  4  and  S  allow  us  to  rewrite  die  first  equation  of  Andreas  as: 


X  t 


Pv(z> 
Py.(T(Z)  ) 


-  1  -  y(r) 


p(r) 

Dv(i)  Pv,(T(z)) 


+  1.004  10* 


P(r)  /5417 
T(2)  \T(2) 


(5.14) 


Widi  these  five  assun^ons  equations  (S.1HS.3)  of  Andreas  are  replaced  by  (5.14),  (5.10)  and 
(S.8),  respectively.  A  closer  inspection  of  the  latter  equations  reveals  diat  (5.8)  and  (5.10)  ccmtain 
only  constants  or  variables  depending  on  the  drc^let  radius  r  and  the  dry  radius  r^  Thus,  ghreo 
die  radii  of  the  parent  and  offspring  categories,  die  parameters  y(r)  and  p(r)  are  fixed. 

S.4  The  source  terms  in  the  SeaChise  model 

Let  us  now  tom  to  the  modeling  of  the  source^sink  term  Sg(r)  for  droplets  widi  radius  r.  Mesti^er 
and  Lefauconnier  [1988]  derived  an  expression  for  die  source-sink  term  S^fir)  for  droplets 
(equatkm  (2.9)  in  chqiter  2): 

Sn(r)  =  6p,  -  (5.15) 


under  dre  conditions  diat  the  total  number  of  dit^lets  in  die  calculation  domain  is  conserved  and 
diat  the  bms  of  tbe  droplets  are  adjacent  hi  die  above  eiqiressitm,  p,  denotes  die  m—s 
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coocentiatioo  of  droplets  widi  ndhis  r.l^l  £a  contrast  to  Ronanlt  etaL  [1991],  adto  sinqdified 
(S.IS)  assuming  diat  die  tenn  rt  (see  (S.14))  dqioids  mily  weakly  on  r,  we  will  not  fmtfaer 
simplify  (S.15). 

Equation  (5.15)  governs  the  exchange  between  dnqdets  of  difEerent  ladiL  Since  the  evaporadan  of 
sea  spay  droplets  widi  ejection  radius  r^i)  is  monitored  via  its  ofEqving  categories,  dnqilets  of 
different  parent  categories  ii  and  ii  are  indqpeodent  Thus,  we  use  (5.15)  to  calculate  die  exdiange 
between  a  parent  category  r^ip  and  its  varioos  offqiring  categories  rg^ii  j),  and  we  iqieat  this 
for  all  n  parent  categories: 


S„(i) 

Sn 


Et"  S„(i.j) 
K..  s„(i) 


(5.16) 


The  total  sum  S,  allows  us  to  calculate  die  source  toms  of  tenqieniture  and  water  vqior,  which 
coiqile  all  indqiendent  parent  categories  (see  (2.3)  and  (2.5)  in  duqiter  2): 


Sy  “  -  Sn 


Sj 


Ly 

Pa  Cp. 


Sa 


(5.17) 


9 


Ite  we  or  ihe  qroAol  for  droplM  cooeeauiiaB  (mUHe)  aqr  bad  10  ODieMioB,  becMue  the  qmboi  p  it  alto  owd  to  doDMe  dewiliM  of 
labitMfW  (ccBtmO.  Howanr,  bDce  dra^  cwcwwibw  la  te  SiaClwt  wodet  an  aapwiaed  ta  (ktbi?],  ihe  qndbol  aacan 
appropUaiB. 
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6  IMPLEMOrrATION  IN  THE  SEACLUSE  C(E»B 


6.1  Calculation  of  radii  of  droplet  cal^ories 

The  radii  of  the  (^E^aiiig  categories  should  be  cfaosea  in  accordance  widi  the  expected  ahiinkiiig 
of  die  droplets,  ududi  is  timtted  to  an  eqnilibriiiin  size  detconined  1^  die  amlnent  relative 
humidity.  Therefore,  Mestayer  [1990]  uses  the  lowest  value  of  die  relative  humidity  in  die 
cakulatkm  dL.uiain  to  infer  the  smallest  radius  of  die  offering  categories  by  the  procedure 
outlined  below. 


At  equilibrium,  f  s  0,  and  (S.14)  reduces  to: 


y 


Pv<z)  -  1  .  RH 
Pv,<T(2)  )  ”  100 


0.6776  |l 


(6.1) 


where  we  have  also  used  (S.IO).  In  die  case  of  equilibrium,  y  and  p  no  longer  dqpend  on  die 
dr<^let  radius  r,  but  only  cm  die  relative  humidity.  Consequently,  the  salinity  and  droplet  density 
are  also  indqiendent  of  r  (see  (5.9)).  The  tn«Tminin  mixing  ratio  in  die  domain  is  given  by: 


Pmx  “ 


1 

0.6776 


I 

100  j 


(6.2) 


Subsequentty,  (5.8)  is  solved  iteradvety  to  yield  p^,  using  (5.9b)  to  calculate  the  droplet  salinity 
S — .  and  wiA  an  initial  value  oipasp^  1>picany,  diis  takes  5  iterations. 


Equation  (5.9)  is  used  once  more  to  infer  a  relatitm  between  die  dry  radius  of  the  particle  and  its 
ejection  radius.  To  diis  end,  it  is  noted  diat  the  original  salinity  and  density  of  die  particle  are 
equal  to  die  salinity  and  density  of  sea  surfeoe  water 


Ps  rlty 

P«w  I«j 


3796 


(63) 


Combining  (6.3)  and  (5.9b),  r^(i)  of  each  parent  category  i  can  be  calculated: 
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Finally,  an  equidisduit  ladiDS  grid  is  built  for  dte  offqnog  categories  (ij)  of  pocnt  categoiy  i: 


Ar(i) 

r(i.j) 


(i)  “  —  /j  \  ^  "  ^mx 

«.«  - 1 - 

r.j(i)  -  <j-l)  Ar(i)  1  s  j  s 


{65) 


Strictly  speaking,  diere  are  only  (iigipl)  ofiFqning  categories.  For  we  have:  i(i>l)  *  or  in 
words,  the  parent  category  is  also  the  first  offspring  category. 


Relation  (63a)  furnishes  die  key  of  the  grid  of  die  parent  categories,  Le.,  their  bin  widths  are 
proportional  to  dwir  nominal  radius.  This  defines  a  non-equidistant  parent  grid:  die  smaller  r^^Ci), 
the  smaller  die  bin  widdi. 

6.2  Calculation  of  r-dqiendent  variables 

As  mentioned  in  die  prevkms  chapter,  D^'(t)  (eq3.11),  p(r)  (eq3.8)  and  p(r)  (eq3.9)  may  be 
calailated  once  the  radii  of  die  drc^tet  categories  have  been  determined.  Unfortunately,  (S.8)  and 
(5.9)  are  coupled  and  have  to  be  conqmted  together  in  an  iterative  process.  Mestayer  [1990]  gives 
die  initial  values  for  diis  calculation: 


Po<r) 


,£± _ 1. 

Pw  (r/r^)*  -  1 

-  Pw  [  1  ♦  Po(i)] 


(6.6) 


lypicalfy,  )t(r)  and  p(r)  converge  after  10  iterations. 

In  princqde,  ii^)  and  p(r)  are  a  function  of  bodi  parent  category  i  and  of^xing  category  j,  diat  is: 


TNOraport 


P«ge 

47 


jl(r)-H(g) 

p(r)-p(i4) 

However,  we  have  choaoi  a  radius  grid  in  ndiidi  11,0  does  not  dqwad  <»  L  hi  that  case  die  ratio 
r(iJ)A'«j(i)  i*  indqiendent  of  L  We  verify  diia  by  rearranging  (6J): 


r(l.j)  .  3^  _  (j-1)  Ar(l)  .  ^  _  (j-1)  (l-W 
x.j(i)  r,3{i)  nott-1 


(6.7) 


Conseqoentfy: 

it<r)«ii<j) 

P(r)»p(j) 

6.3  T||y  and  Vf  for  droplet  categories 

Values  for  die  flight  tune  Tg^CiJ)  and  final  fidl  velocify  V|(iJ)  are  retrieved  from  a  database,  hi 
general,  die  database  will  not  have  entries  for  the  radii  r(ij)  and  Unear  interpolation  of  die  nearest 
entries  r^  and  r.  is  dien  used.  When  die  radios  r(ij)  is  thu  die  smallest  radius  rgig  in  die 

database,  die  final  fall  velocity  is  ^^oximated  by  die  Stdres  Velocity: 


Vj(i,j) 


(6.8) 


This  ifiproximation  is  onfy  correct  for  r  <  SO  pm  [Rooauh  etaL,  1991],  and  SeaChise  issoes  a 
warning  when  r  >  50  pm.  When  r(ij)  <  r^,  the  fli^t  tune  is  extrqiolated  from  die  first  two 
entries  ri  and  r2  in  die  database: 
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"5r  r,  -  ri 


(6.9) 


6.4  Initial  and  boondaiy  cmiditiofu  tor  SeaOuae 

The  initial  conditions  for  die  ofEqving  categories  should  be  identical  to  the  initial  conditions  for 
the  parent  categories,  hi  diis  manner,  die  model  run  starts  widi  small  (no)  gradients  dpjdt  (see 
(S.IS)).  Therefore,  we  use  die  initial  condition  previously  given  by  Rouault  etaL  [1991]  (see  also 
chapter!): 

p,(i.j,z)-0  Vi.j.,  (6.10) 

We  have  not  changed  the  boundary  conditions  at  die  top  and  die  bottom  of  the  calculatioo 
domain,  diat  is,  we  use  die  equatkms  giveo  in  charter  2  for  bodi  the  parent  and  die  oCfqxing 
droplets: 


Px  (i»  j  *0 


(6.11) 


P,(i»l/Z.in) 
P*  (i»  j 


jr^s„(i)dz  ♦  p.3(i>  v^(i) 
vTTTTT) 


/^S,(i,j)dz 

"  *^Ve(i,j)  2  i  j  « 


(6.12) 


where  p^i)  denotes  die  concentration  oS  ejected  dropleu  with  radhu  r^^)  and  we  have  used  diat 
p^  *  0  for  offiqiring  catteries.  The  boundary  condition  ditt  die  gradient  in  pr  must  be  zero  at  die 
top  of  die  domain  is  not  nupoeed,  but  is  dieefcad  at  die  end  of  die  model  tun  (see  below). 
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6J  ThemodclpifOTewwCiMdCainSoClme 

Once  we  me  eqiesnooi  given  previomfy  by  [Rananlt  etiL,  1991]  to  introdoce  the 
coimterdiffmion  and  die  meooacoiac  flnx  (aee  eqmtiooi  (2.7)  and  (2.8)),  atfafMed  for  patent  and 
oCE^niig  categofiea: 


CiTjiyd,!)  ^  -  P,(i.z)  ] 

1  (6.13) 

CiTfiyd,  j)  ^  1  2  s  j  s  no« 


SC^Ci,  j)  *  SCt  (  1  ♦  C,  y*----r  1  1  (6.14) 

1.56  U, 


edwte  pP  denotes  die  concentratkm  profile  of  parent  cat^ory  i  in  non-tnrtwlent,  non-evqxnative 
conditions  (see  ch^iter  4  and  below).  Obvioosty,  p,®  «  0  for  off^xing  cat^oties  since  diey  onty 
get  populated  by  evqioFstion. 

6.6  The  droplet  profiles  in  a  non-tuibalent.non'evqiocativeatinoqiliere 

The  vertical  profiles  of  droplet  concemratko  in  o«i>tiirbalent  and  non-evqNxative  conditioos 
are  used  in  die  modeHiiig  of  die  macroscopic  flnx  of  die  droplets  dne  to  ejectioo  and  gravitational 
fidl  (see  (2.8)  in  dupter  2).  Notmatised  profiles  (^fi,z)  have  been  calculated  previously  (aee  (4.5) 
and  are  available  in  a  database.  The  nmidiers  CCM)  are  converted  into  coocemration  Pr'^i.z)  in 
by  mult^lkation  fvith  die  density  p^)  of  droplets  ejected  at  die  wave  snftce: 
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so 


Px(l,z)  ■  C(i,*)  p,j(i)  (6.15) 

The  density  p^)  is  found  from  die  lelsdoo  between  the  scwrce  foncdon  dF/dt  in  [arV>|im'>]  and 
die  flu  ejected  drofilets  p^)V^(i): 


.  Ar(i)  .  A,p^>j(i) 


(6.16) 


V^(i) 


edicfe  Ar(i)  is  die  size  bin  of  parent  categcoy  L  ^  a  first  guess,  the  source  function  dF/3r  is  taken 
from  Andreas  [1992],  who  uses  data  from  Miller  [1987]  for  small  droplets  (0.8  pm  ^  r|o  ^  IS 
pm): 


Bo  ♦  Bjdog  r,o)  ♦  B,(log  r,o)* 
♦  B3(log  r,o)’  +  Bo  (log  r,o)* 


(6.17) 


and  data  from  Wu  etaL  [1984]  for  larger  droplets  Otiq  >  IS  pm); 

15  <  Zgg  <  37 .5  pm 

37.5  i  Tgo  S  100  pm  (6.18) 

z,o  k  100  pm 

In  these  source  fhnetioos,  rgg  denotes  die  radiiis  of  the  dio|det  at  an  ambioit  relative  humidity  ctf 
80%.  To  convert  rgg  into  r,  Andreas  [1989, 1992]  gives  the  following  conversion  formula's: 
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r,o  -  0.518  r®  ”‘ 


dF 


dlao  dF 


0.506  t-® 


(6.19) 


The  ooefiSdeots  B0-B4  in  (6.17)  end  CyC^  in  (6.18)  »  a  function  of  Uk).  Andreas  [1992] 
tabulates  values  for  dteae  coefficients  in  die  wind  qieed  range  6-18  in/s.  By  mteqwlatkm,  Sea- 
Qnse  determines  die  value  ci  dPidt  tat  the  value  of  U]o  qiedfied  for  die  model  lun.  Since 
extr^lation  is  not  (yet)  suppotted,  the  range  of  values  td  Uiq  in  SeaGuse  is  restricted  to  6-18 
m/s. 

The  source  function  given  above  is  by  no  means  a  consensus  source  function.  The  establishment 
of  such  a  source  function  is  rather  difficult,  because  ejqxarimental  data  are  i^arse  and  hard  to 
accompli^  An  elaborate  discussion  of  diis  subject  would  be  outside  die  scope  of  diis  iqwct  (see 
Andreas  [1992],  and  references  dierein).  However,  the  validity  of  Andreas'  source  function  is 
questionable  [Katsaros  and  De  Leeuw,  1993].  The  two  major  objectkHis  are; 

a)  The  dataset  of  Miller  is  small  and  shows  a  large  amount  of  scatter  in  the  regioo  105rgo5  IS 
pnL 

b)  The  dataset  of  Wo  is  based  on  data  obtained  in  low-wind  qpeed  cmiditioos  (Uio<  8  m/s)  and 
therefore  not  rqvesentative  for  spume  droplets. 


6.7  Calculation  of  the  saturated  water  vyir  density 

In  order  to  calculate  die  saturated  water  vapor  density  p„(T)  in  (5.14),  we  use  Bock's  formula 
[1981]  to  compute  die  saturated  water  viqwr  ptessure  e^: 

/  17. SOS  a  \ 

•aat^T)  ■  (1.0007  ♦  3.46  10*  P)  6.1121  e'  (6.20) 

where  0  ■>  T  -  273.15,  and  die  standard  pteswite  ^1013  mbar  is  used  to  obtain  e,JT)  in  mbar. 
Subsequently  we  die  ideal  gas  law,  as  suggested  by  Andreas  [1989]  to  conqime  6ie 
saturated  water  vqwr  density: 
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Pv,(T) 


100  e^t(T) 

R  T 


(6^1) 


widi  T  in  K  to  obtain  (\,(T)  in  [kg/ioa^]. 

6.8  The  vettkal  grid 

The  vertical  grid  in  the  SeaCltiae  model  hat  to  extend  from  dte  wave  troa^  to  die  top  of  die 
boundaiy  layer.  As  mentioned  ptevioosly,  we  have  divided  the  vertical  axis  into  two  xegitms:  a 
wave  region  below  the  wave  crests  and  an  atmoqdieric  region  above  die  wave  crests.  The  origin 
of  die  vertical  axis  is  located  at  mean  sea  leveL 

The  grid  in  the  wave  region  starts  at  a  hei^t  Zq  above  the  wave  trough  (located  at  z  >  -0.46  Hg). 
Our  previous  woik  on  the  CLUSE  model  for  fresh  water  and  die  first  tests  of  die  SeaChise  model 
have  shown  that  it  is  necessary  to  start  widi  a  small  grid  spacing.  We  found  that  the  spacing 
should  be  of  die  order  of  the  value  of  the  eddy  diffiisivity,  or  die  numerical  solver  wiU  not  work. 
Since  the  eddy  diffiisivity  near  die  wave  trough  is  radier  small  and  the  nuinber  of  grid  points  is 
limited  for  conqmtational  reasons,  a  logarithmic'  grid  [Rouault  et  al.,  1991]  of  SO  grid  points  is 
introduced: 

IZ[1]  »  -0.46  Ho  +  2o 

,  (6.22) 

z[il  *  z[i-l]  +  A2[1]  (2  s  i  s  50) 

where  Az[l]  is  typically  in  die  order  1(H  m  and  die  irregularity  parameter  x  is  diosen  such  that 
Az[30]  is  in  die  order  of  1-S  centimeters. 

The  remainder  of  die  wave  region  is  coveted  by  a  linear  grid,  widi  a  grid  qmcing  equal  to  Az[S0]. 
In  this  manner  we  ensure  diat  die  regkm  below  the  crests,  uhere  die  productkii  of  sea  Mpny 
aerosol  takes  place,  is  well  coveted  by  die  grkL  Continuing  die  non-linear  grid  would  lead  to  an 
unacceptable  large  grid  ^pacing  near  die  crests.  Thns: 


z[l]  -  z[i-l]  *  Az[50] 


(51  S  1  S  N,^) 


(6.23) 
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irfiere  is  die  total  namber  of  grid  pomts  in  the  wave  regk».  To  esirare  that  z[NavJ 
pnqiriAty  eucdy  with  die  wave  oeat.  the  ittetval  between  z(NaM«-ll  ood  is  somewhat 

enlarged  orreduced. 


The  atmoqdieric  legkn  of  die  grid  fron  die  wave  crests  (at  z  ■  O^SHq)  iq>  to  die  top  of 

die  cateulation  domain  (at  z  »  haJ.  Since  die  exchange  of  vapor  of  die  freshly  generated  sea  qxay 
aerosols  mainfy  takes  place  in  die  lower  regions  of  die  atmo^dieie,  it  seems  appn^ciate  to  use  a 
grid  in  die  atmoqdieric  regkm  which  is  more  dense  near  die  wave  crests: 


2[il  «  z[i>l]  +  X*'*  Az  [Naav.+1] 


A2lNa,a.+l] 


(htot-0>53Ho)  (1-x) 
i-x*'*- 


(6^) 


uriiere  die  number  of  grid  points  N,aa  ^  •djusted  to  obtain  an  initial  spacing  Az[Naa,e'*’l]  that  is 
close  to  AzlSO}.  Typical  values  for  the  irregularity  parameter  %  sre  1.04-1.08,  leading  to  a  grid 
spacing  of  1-5  m  at  z  100  m.  The  total  number  of  grid  points  (N«ne  ‘t’  is  then  in  die  order 
of  200. 

6.9  Checks  during  die  SeaQuae  model  run 

During  die  model  run,  duee  checks  widi  regard  to  die  calculation  of  die  source  terms  S,,  S,  and 
Sx  are  made: 


>1:  f(ij)>0  (8eeeq3.14) 

This  means  diat  ditiJVdt  >  0  and  that  die  droplet  tends  to  condensate,  lii  odier  wtnds,  die 
SeaCbse  model  moves  away  from  the  correct  solution.  However,  in  die  end  die  model 
should  still  be  able  to  arrive  at  the  correct  solution,  provided  diat  die  timestqp  of  die  model  is 
sufiBciently  smaU.  Therefrne,  die  case  f  >  0  is  noted,  bat  no  frirdier  action  is  taken. 


>2:  p,(ij)<0  (seeeq.5.15) 


AMwugh  negative  oonoentrationa  for  any  droplet  categoey  aie  physically  inq)08sd>ie,  we 
catmoc  just  set  p,(y)  to  0.  hi  diat  case  we  would  violate  die  model  ocmatraint  diet  the  number 
of  droplets  is  conserved.  At  dK  end  of  die  model  ran,  all  p/ij)  duNild  be  positive  (again). 
Thetefose,  die  case  pj  <  0  is  noted,  but  DO  tedier  action  is  taken. 

►3!  (s  0.  4  -  p.  ^  )  >  0 

The  mmerical  mediod  (see  chapter  2)  used  in  die  modd  has  the  constramt  that  diis  part  of 
(S.15)  must  be  negative.  Therofore,  in  die  case  that  the  above  factor  is  positive,  it  is  noted 
and  die  factor  is  reset  to  zero. 

6.10  Checks  at  die  end  of  die  SeaGnse  model  run 

At  die  end  of  the  model  run,  duee  inyottant  checks  idaied  to  die  validity  of  die  final  results  are 

made: 

>1:  The  value  RH_i_  of  the  relative  himaiditv  in  the  domain  is  determined.  If  is  much 
lower  than  die  value  previously  determined  to  set  up  die  offqxing  radius  grid  (see  $6.1),  die 
radius  grid  may  not  extend  to  sufficiently  small  radii  and  droplets  may  not  have  been  able  to 
evifxirate  to  their  equilibriiim  sizes,  hi  such  a  case,  the  calculatioD  should  be  restarted  widi 
an  extended  offspring  radius  grid. 

>2:  The  boundary  conditioa  dp,/9z  »  0  at  die  top  of  the  calculation  domain  (see  chi^iter  2)  is  not 
inq  oaed  during  die  model  run.  Tbeiefore,  it  must  be  checked  at  die  end  of  die  model  run  to 
ensure  that  the  number  of  droplets  in  die  domain  has  been  conserved. 

>3:  The  cwvergence  of  die  model  solution  is  verified  for  die  profiles  of  water  vqior, 
temperature  and  droplets.  Coovetgeoce  is  checked  indqieodratfy  for  water  vt^xv, 
tenqieratnre  and  droplet  otmcentratioii. 

6.11  I^esmitatioottf  final  results 

One  of  our  aims  is  to  compare  the  resuits  of  the  SeaChise  model  with  eiqieritiiental  data  obtained 

widi  (optical)  partide  coanters  or  Rototods.  However,  it  is  not  posrible  to  directly  conqiare  die 

model  results  and  die  data,  becwse  die  experimental  data  is  presented  in  adjacent  size  Ians 
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R|(k)  tViARfCk)  and  die  model  lesultt  in  overiapping  parent  and  offspring  categories.  In 
yAHtirm,  die  dropkt  concentnakos  in  die  SeaQnse  model  are  expressed  in  [kgAn^],  ^idmeas 
experimesital  data  is  often  presented  in  [pn^/cm^. 


Hwrefore,  the  parent  and  o£Eqiriiig  categories  r(ij)  ±  Vidr(i)  mast  be  nupped  onm  dw  radius  grid 
R«(k)  ±  ViARfCIt),  <H>ce  die  model  nm  is  completed.  The  nupping  is  done  by  transferring  a 
prapoitioiial  part  of  p^ij)  to  Rf(k)  wbeoever  the  size  bin  of  droplet  category  r(ij)  (partially) 
overlaps  with  a  size  bin  Rf(k).  After  diis.  we  convert  the  p,  values  (in  kgAn^)  into  dV/dr  values  (in 
pm2/cm3)  by: 


dV 

■S’ 


Pw  Ar 


(6.25) 


where  Ar  is  in  pm. 


TNOraport 


Page 

56 


7  DISCUSSION  AND  AIMS  OP  FUTURE  RESEARCH 

The  previous  clusters  have  outlined  the  SeaQuse  model,  which  extends  the  CLUSE  model  to 
'(^>01  ocean'  conditions.  Several  new  hypotheses  had  to  be  made  in  order  to  describe  the  air  flow 
ova  the  waves  (chqiter  3): 

1.  The  waves  are  modelled  by  a  2-dimensional  Stokes  wave 

2.  The  waves  are  fully  developed 

3.  Momentum  is  only  transferred  to  the  water  (no  form  drag) 

4.  The  atmos{Aeric  stratification  is  neutral 

An  additional  hypothesis  was  introduced  to  describe  the  movement  of  the  droplets  in  the  air  flow 
(chipter  4): 

5.  The  droplet  flight  times  in  non-turbulent  and  non-evaporative  conditions  are  finite 

FurthemKne,  approximations  were  made  in  the  equations  that  govern  the  eviporadon  of  salt  water 
droplets  (chapter  S).  These  limit  the  applicability  of  the  SeaCluse  model  to: 

r>lpm  0ST^30*C  RH^75% 

Whether  or  not  the  new  hypotheses  are  accq>table  cannot  be  discussed  here,  because  the 
SeaCluse  model  has  not  yet  beat  tested  and  validated  with  experimental  data  from  the  HEXOS 
program  [De  Leeuw,  1990a].  Testing,  and  subsequent  conparison  of  the  model  results  with 
experimental  data  will  be  the  primary  objective  in  the  continuation  of  this  research. 

The  outcome  of  these  tests  determines  whether  our  model  for  the  air  flow  over  the  waves  needs  to 
be  further  elaborated.  Also,  it  may  be  necessary  to  modify  or  to  replace  the  aerosol  droplet  source 
function  (see  the  discussion  in  chipter  1  and  6.6).  New  and  more  detailed  experimmtal  results 
on  the  number  of  jet  droplets  released  from  bursting  bubbles,  their  ejectitm  heights  and  their  velo¬ 
cities  will  be  available  soon  [Spel,  19928  and  1992b]  and  may  be  included  in  the  SeaCluse 
model. 
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The  SeaChue  model  only  takes  jet  droplets  into  account,  neglecting  film  (boplets  that  oe  also 
produced  fiom  bursting  bubbles.  In  fact,  die  larger  bubUes  produce  no  jet  drops  at  all,  ooty  film 
droplets  [Blanchard,  1983].  Even  so,  die  relative  cootributimi  of  film  drops  is  small  [Wu,  1989], 
because  die  film  droplets  are  very  »miiH  and  are  ejected  with  small  vertical  veloctties.  Therefore, 
we  feel  diat  neglecting  die  film  droplets  is  justified. 

hi  addition  to  the  bubble-mediated  production,  qmme  droplets  may  be  direcdy  produced  wave 
tearing.  This  mechanism  comes  into  pls^  at  wind  qieeds  exceeding  9  m/s  [Monahan  etal.,  1986]. 
The  question  as  to  what  extent  die  spaac  droplets  contribute  to  die  total  dn^let  flux  has  not  yet 
been  definitively  answered.  Wu  [1990]  estimates  a  contribution  of  33%  at  wind  qieeds  of  13  m/s, 
based  <ni  experimental  by  De  Leeuw  [1986].  On  the  other  hand,  more  recent  experimental 
evidence  by  De  Leeuw  [1990b]  shows  diat  i^iume  dn^lets  are  unlikely  to  contribute  much  to  die 
total  aerosol  concentrations  at  tins  wind  speed. 

The  SeaCluse  model  may  be  extended  widi  a  module  that  describes  the  production  of  qiume 
droplets  at  die  wave  crests,  hi  diis  manner,  the  effect  of  spume  drops  <»  the  humidity  profile  and 
the  heat  fluxes  can  be  studied.  However,  it  should  be  emphasized  diat  only  qualitative  results  are 
to  be  eiqiected  in  die  absence  of  a  reliable  aerosol  source  fimction. 

Our  model  of  die  air  flow  over  the  waves  leads  essentially  to  a  wave-rotor  between  die  wave 
crests  (conqiare  figures  12  and  3.7).  The  wave  rotor  model  has  bem  proposed  De  Leeuw 
[1986]  to  explain  die  observed  mininaa  and  maxima  in  die  particle  concentration  profiles  (10  <  D 
<  100  pm)  between  0.2  and  20  m  AMSL  [De  Leeuw,  1986,  1987].  With  our  wave-rotor  modeL 
we  conjecture  diat  the  largest  dn^let  crmcentrations  will  exist  in  die  lee  side  of  die  waves.  This 
conclusion  is  based  on  die  extended  droplet  suspension  times  at  diat  locaticm  and  die  droplet 
trajectories  (see  figures  4.1  and  4.2). 

This  is  in  accordance  widi  laboratory  eiqieriments  by  Koga  and  Toba  [1981],  who  also  observed 
larger  concentrations  at  die  lee  side  of  the  waves,  resulting  fiom  drt^lets  produced  at  die 
downward  shqie  of  the  same  wave  trough.  On  the  odier  hand,  De  Leeuw  [1989c]  indicated  diat  in 
die  field  die  production  may  be  higher  on  die  windward  side  of  die  crest,  because  the  uhitecuqi  of 
a  trailing  wave  can  be  at  the  windward  side  of  die  crest  This  couhi  be  an  inqietus  to  include  a 
module  in  SeaOuse  that  describes  dw  wave  breaking. 
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The  SeaCbiae  model  yields  infonnarion  oo  the  vettkal  ttuctme  of  the  atmoqthae  (water  vqior, 
droplet  dtstnbotMm,  tetnperatnre),  from  wfaidi  myortanf  parameters  sodi  as  die  refractive  index 
and  die  extinction  may  be  infetred.  Therefore,  die  model  could  be  used  to  assess  the  inaxiimim 
effective  range  of  electro-optical  sensors.  One  of  the  requirements  for  an  operational  range 
predictor  is  diat  die  input  parameters  can  be  easify  obtained.  In  diis  reqwct,  the  SeaCluse  model  is 
nearly  perfect:  its  hqiot  consists  of  onfy  four  parameters,  diat  can  be  obtained  by  standard  instru¬ 
mentation:  U|Q  KHjq  Tjq  I'tKf* 

In  addition,  die  SeaQuse  model  uses  the  drqilet  source  function  givoi  by  Andreas  [1992].  This 
function  (see  $6.6)  is  a  function  of  U|o  only. 

Theoretically,  die  smallest  droplets  that  can  be  handled  the  SeaCluse  model  have  a  radius  of  1 
pm.  However,  computational  demands  become  radier  hi^  udim  die  droplet  radius  is  in  the  order 
of  a  few  microns.  This  is  a  severe  drawback  for  die  iq^licatimi  of  the  model  in  range  prediction, 
since  droplets  widi  radii  smaller  dian  10  pm  determine  the  visibility  in  die  visible  wavelength 
region  and  contribute  appreciably  to  die  visibility  in  the  IR. 

In  order  to  overcome  diis  drawback,  die  CLUSA  model  is  being  developed  at  the  Ecole  Centrale 
de  Nantes  [Mestayer  etaL,  1991].  CLUSA  describes  die  dupersion  of  a  distributimi  of  aerosols 
from  0.0S  to  2.S  pm  radius  in  the  marine  atmospheric  surface  layer.  In  its  present  state,  CLUSA 
starts  from  an  initial  profile  of  dry  aerosol  dutt  is  transformed  to  marine  conditions  by  calculating 
die  aerosol  growth  by  water  vapor  condensation.  Subsequently,  die  effective  dqmsition  velocity 
of  the  aerosols  is  calculated  as  a  functirm  of  dieir  ori^nal  size  and  chemical  comporiitimi,  of  die 
bygro-dietmodynamic  atmospheric  condition,  and  of  die  sea  state.  The  calculation  includes 
gravitational  settling.  Brownian  difiusun  and  the  presence  of  the  diffusive  sub-htyer  close  to  the 
sea  surface. 

bn  the  future,  merging  CLUSA  and  SeaChiae,  die  model  should  include  die  interactions  of  the 
smaller  aerosols  widi  sea  spray  dxopkiu  by  coalescence  and  widi  die  water  vapor  field.  The 
condiined  modd  would  dien  govern  dm  complete  size  range  of  aerosols  diat  determine  the 
propagation  properties  of  die  atmoqdiere. 

The  SeaCluse  model  onty  deals  widi  fredily  produced  sea-qaay  aerosols.  The  total  aerosol 
concentration  in  die  atmot;diere  results  from  a  balance  between  local  production  and  removal  and 
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Ae  advectum  of  aetoaol  produced  eltewbere.  &  it  oonjectmed  diet  die  advected  aerostd  is  in 
effect  a  ’backgrmmd*  that  adds  to  the  profile  doe  to  the  locally  prodnced  aeroaoL  hi  view  of  the 
long  mixing  times,  die  latter  constitutes  onfy  a  small  firactioo  of  die  total  aerosol  concentration. 

The  background  aerosol  consists  of  a  mixture  of  man-made,  biological  and  natnral  aerosfd  of 
continental  and  marine  txigin,  containing  a  large  variety  of  hygroscopic  and  non-liygroscc^  tpt- 
cies  ndiich  strongly  influmce  die  dependence  on  relative  humidity.  Ftirdiennore,  aerosol  dropl^ 
are  not  inert  and  their  cootyositioo  can  diange  through  many  chemical  reactions  iriiich  take  place 
in  the  atmosphere. 

Van  Eijk  and  De  Leeuw  [1992]  recently  addressed  the  inqiattance  of  taking  die  background 
aerosol  of  continental  origin  into  accotmt  uhen  predicting  propagation  parameters  (e.g., 
atmospheric  extinction)  for  a  coastal  region,  .^iparently,  well-establidied  aerosol  models  such  as 
the  Navy  Aerosol  Model  [Gatfaman,  1983]  fiul  for  ^tinted  inner  seas'  like  the  North  and  the 
Mediterranean  [Tanguy  etaL,  1991]. 

hi  conclusion,  a  model  that  aims  to  provide  an  adequate  descttytion  of  the  propagation  properties 
of  the  atmosidiere  cannot  be  confined  to  die  fieshty  prodnced  sea-salt  aerosols.  In  addition, 
modules  diat  handle  the  continuously  changing  cbonical  conqiosition  and  die  diqiersion  of  die 
'background'  aerosols  are  needed.  Onty  dien  die  model  can  truly  govern  the  production,  advecticm 
and  dqiosition  of  the  aoosols. 
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NOMENCLATURE 


A.1 


CLATURE 

Symbol 

Unit 

Ow 

• 

water  activity  an  aqueous  sohition 

c 

mol/1 

molar  concentration  of  an  aqueous  solution 

V 

J/kgK 

qtecific  heat  of  air  at  constant  pressure 

CjJN 

Drag  coe^ient  (noitral  stratification) 

C, 

m/s 

phase  velocity  of  die  waves 

D 

m 

D, 

m^/s 

molecular  diffiisivity  of  water  vapor  in  air 

d; 

m2/s 

modified  Dy 

g 

m/s2 

acceleration  of  gravity 

htot 

m 

height  of  the  calculation  domain 

Hcj 

m 

jet  droplet  ejection  height 

Ho 

m 

wave  amplitude 

i 

index  of  pareid  categories 

j 

index  of  offspring  categories 

k 

m-* 

propagation  constant 

K 

W/mK 

modified  diennal  conductivity  of  air 

Lv 

J/kg 

latmt  heat  of  vqxirizaticMi  of  water 

m 

kg 

mass  of  die  dn^let 

m, 

kg 

mass  of  NaCi  in  die  droplet 

niw 

kg 

mass  of  pure  water  in  die  droplet 

M. 

kg 

molecular  weight  of  air 

M, 

kg 

molecular  weight  of  NaQ 

Mw 

kg 

molecular  weight  of  water 

n 

number  of  parent  categories 

noff 

number  of  offiqning  categories 

Nx 

number  of  hcoizontal  grid  points 

Nz 

number  of  vertical  grid  points 

ft. 

Ihandtl  number 

r 

m 

radius  of  droplet 

m 

dry  tadhu  of  die  droplet 

m 

radius  of  parent  category  (ejection  radius) 

'off 

m 

radius  of  offiqning  category 

R 

JAnoUC 

universal  gas  cmistant 

Re 

Reynolds  number 

RH 

% 

relative  humidity 
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NOMENCLATURE  A.2 


RH«f 

% 

RH,o 

% 

S 

s„ 

S<|* 

Sv 

s,* 

mM/ 

mol3/2 

Sc, 

T 

K 

Tfly 

s 

Tjnif 

•c 

T,o 

•c 

T. 

•c 

U 

m/s 

Ujo 

m/s 

m/s 

«r 

m/s 

U« 

m/s 

V. 

afi/mol 

V,0 

m^/mol 

Vf 

m/s 

m/s 

W 

m/s 

Wd 

m/s 

Wr 

m/s 

t 

s 

X 

z 

m 

zo 

m 

ZOL 

m 

zo. 

m 

Oc 

Or 

Y 

8p 

m 

At 

m 

Aw 

m 

relative  humidity  at  die  surface 

relative  humidity  at  z  »  10  m 

salinity  of  the  droplet 

soutce-smk  term  for  droplets 

salinity  of  sea  surface  water 

source  tenn  of  temperature 

source  term  of  water  vapor 

function  of  T  in  finding  die  apparent  molar 

volume  of  an  aqueous  solution 

Schmidt  number 

temperature 

flight  or  suqiension  time  of  drc^let 
sea  surface  temperature 
temperature  at  z  »  10m 
tenqimture  scaling  factor 
horizontal  wind  ^leed 
horizontal  wind  qieed  at  z  »  10  m 
horizontal  speed  of  droplet 
horizontal  speed  of  droplet,  relative  to  air 
friction  velocity 

rqiparent  molal  volume  of  an  aqueous  solution 

V,  at  infinite  dilution 

final  fall  velocity  of  the  droplet 

droplet  ejection  velocity 

vertical  wind  speed 

vertical  speed  of  droplet 

vertical  q>eed  of  droplet,  relative  to  air 

time 

horizontal  comdinate 

vertical  coordinate,  z  =  0  at  mean  sea  level 

roughness  length  calculated  with  Chamock  <»' 

Kitaigorodsld  nxrtel 

Zq  for  local  profiles  of  U,  p,  and  T 

Zq  for  the  atmospheric  profile  U(z) 

empirical  constant  used  in  corrqniting  D,' 

empirical  constant  used  in  conqwdng  Ic,' 

total  number  of  ions  into  which  a  salt  molecule  in  die 

droplet  dissociates 

step  size  in  droplet  trajectory  calculation 
empirical  length  scale  used  in  conqmting  k,' 
empirical  length  scale  used  in  computing  D,' 
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A.3 


At 

m 

size  bin  of  droplet  category 

At 

8 

time  8tq>  in  droplet  trajectoiy  raii«^ii«rif»n 

tl 

m 

hei^  of  d»  wave  surface 

K 

Von  Kannan  constant 

X 

m 

length  of  the  wave 

mass  mixing  ratio  of  NaCl  and  H2O  in  the  droplet 

V 

m?/s 

kinematic  viscosity  of  air 

P 

kgfw? 

density  of  droplet 

P. 

kg/m3 

density  of  dry  air 

Pd" 

kgAn3 

concentratum  of  droplets  of  diameter  D  in  non-ev^)orative  and 
nm-turbuknt  conditions 

Pq 

kg/m3 

concentratirai  of  droplets  ejected  from  the  surface 

Pr 

kgAn3 

concentrati(»i  of  droplets  with  radius  r 

P. 

kgAn3 

density  of  NaCl 

P«W 

kg/lm^ 

density  of  sea  surface  water 

Pv 

kg/m3 

water  vapor  doisity 

Pvs 

kgAn^ 

saturated  water  ygpoi  density 

Py* 

kg^3 

scaling  factor  for  humidi^  profile 

Pw 

kgAn3 

density  of  water 

J/m2 

surface  tension  of  a  flat  surface  with  tfy*  aimw 
salinity  and  tenqierature  as  the  droplet 

Ot 

m2/s 

eddy  diBusivity 

♦ 

practical  osmotic  coefficient  of  the  droplet 

atmospheric  stability  function 

X 

8-> 

irregularity  of  votical  grid 

CQ 

wave  period 
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THEWAVEGRID  B.l 


THE  WAVE  GRID 


When  a  variable  is  to  be  averaged  over  dte  wave  smface,  a  grid  is  needed  diat  is  equidistantly 
spaced  over  the  wave  surface.  This  is  not  die  same  as  a  grid  equidistantly  qiaced  in  z  <n  x, 
because  the  wave  surface  has  to  be  follov«d.  When  we  denote  die  coordinate  along  the  wave 
surface  by  s,  the  trajectory  along  die  wave  surface  may  be  calculated  by  a  triangular  mediod: 


Az 

Ax 


tan  a 


cos  a 


Ax 

As 


(B.l) 


Thus,  for  sufficiently  small  triangles  (say,  1000  over  a  wave  period),  we  may  calculate  As  widi 
(B.l),  using  die  derivative  of  the  Stokes  wave  (see  equadcm  3.3): 

*  -TT  sin  (kx-«t)  -  sin  (2  (kx-wt) )  (B.2) 

ax  20  400 

The  total  path  length  S  along  the  wave  surface  is  equal  to: 

S^ZAs 

If  we  neglect  the  curvature  of  the  wave  by  taking  an  x-grid,  the  path  length  X  is  given  by: 


XsZAx 


We  calculated  bodi  X  and  S  and  it  turned  out  that  s  is  only  0.6%  longer  than  the  horizontal  path 
lengdi  X.  Therefore,  we  make  only  a  minor  eizor  if  we  replace  die  wave  surface  grid  (along  s) 
a  horizontal  grid  (along  x).  The  {diysical  reason  for  this  is  diat  die  lengdi  of  die  Stdces  wave  is  20 
times  its  anqilitude,  i.e.,  the  slope  of  die  wave  surface  is  radier  low. 
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EDDY  DIFFUSIVITY  IN  THE  WAVE  ROTOR  APPROXIMATION 


When  we  assume  diat  die  wave  rotor  plays  an  dominant  role  below  the  wave  crests,  die  eddy 
diffiisivi^  in  diat  region  is  gtvoi  by: 


(a^(z,t)  -  Oj 
a^(z,t)  «  0 


(cxest)  z  >  1)  (t) 
2  <  T)  (t) 


where  11(0  is  die  instantaneous  wave  height  at  time  L  To  calculate  o,  below  die  wave  crests,  we 
need  to  know  t4.(z),  i.e.,  the  amount  of  time  diat  height  z  is  above  die  wave  surface  (relative  to  die 
wave  period  Tq). 

We  may  obtain  t^(z)  firom  the  expression  for  the  Stokes  wave  (eq.3.3).  To  this  end,  we  set  x  =  0, 
v^ch  yields: 


11  (t)  = 


-  cos(-«t)  +  cos(-2«t) 

«  80 


Subsequmdy,  this  equaticm  is  inverted  for  0  ^  t  ^  liTg.  Using  die  goniometric  rule  cos  2p»2 
cos^  P  - 1  and  the  abc-tule  for  quadratic  eiqnessions,  diis  yields: 


t*(z) 


± 


azccos 


*  IJ 
20 


The  eddy  diffusiviQr  below  die  wave  crests  is  dien  given  by: 


Ot(z)  »  t*(z)  K  Ot  (crest) 


z  <  fHo 


An  exanqde  of  diis  calculadon  is  shown  in  figure  C.l.  Below  the  wave  crests,  die  ed^  difiiisivity 
follows  the  inverse  of  die  Stokes  wave,  whereas  above  die  waves  die  linear  reladcm  (eq.3.10) 
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EDDY  DIFFUSIVITY  IN  THE  WAVE  ROTOR  APPROXIMATION  C2 


^iplies.  There  is  a  «niMll  diaconthniity  at  the  wave  oeat  caused  by  rounding  emns.  I.e.,  at  die 
wave  atgt  Ot  as  calculated  by  the  linear  function  is  not  exactfy  equal  to  o,  as  calculated  by  the 
inverse  function.  Typically,  die  difference  amounts  to  03%. 


OX)  ai  02  OJ  0.4 


Rguie  Cl:  Eddy  diffesivity  calculated  with  the  wave  rotor  appraximation 
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CHARNOCK  RELATION  AND  KITAIGORODSKI  MODEL  D.1 


CHARNOCK  RELATION  AND  KITAIGQRODSKI  MODEL 

The  SeaChue  model  caknlates  the  foughness  kngih  Zq  ^  CSumock's  leladoa  ox  die 
Kitaigorodksi  model  for  veiificatiao  puipoaes  (see  §3.4). 

Chamock's  leladoii  [Chamock,  1955]  is  giveo  by: 


where  a^  »  0.011  for  open  ocean  [Smith,  1988]  and  a^  »  0.019  for  coastal  areas  (average  value 
firom  [Garatt,  1977],  [Wu,  1980]  and  [Geemeart «  aL,  1986]). 

For  open  ocean  ctmditions,  we  may  also  calculate  the  rou^mess  length  by  die  Kitaigorodski 
model  [Kitaigorodski  et  aL,  1973],  ghrea  by: 


Zo  »  Ak  [  S(«)  e3q?(-2itg/«u,)  dw  j"'*  p2) 

«diexe  «  0.028  [Geeniaert,  1988]  and  S(Q>)  denotes  the  wave  qiectnim.  In  our  modeL  S(o>)  is 
taken  to  be  a  Phil^  wave  spectrum: 

S(«)  ■  p  g*  «■*  «  >  «Q  (D-3) 

where  (Oq  denotes  die  frequency  of  the  Sudces  wave  (see  equatkm  3.3)  and  p  is  given  by  Geemeart 
et9L,[1986]: 

P  »  0.005  ♦  0.002Q  ♦  1 . 5  (g/WoU.)  (D-^) 

for  tqpen  ocean  conditums.  Note  that  this  relation  for  p  was  established  for  wave  qiectra  measured 
during  the  MARIEN  oqieximent  for  unlimited  fetch  conditioos  but  for  water  dqilhs  of  onty  15 
meter.  For  qiplicataon  to  coastal  seas,  dus  model  should  be  inqaoved  to  take  into  account  the 
limited  fetch  in  die.case  of  wind  from  die  shore. 
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CHARNOCK  RELATION  AND  KITAlGCHtODSKl  MODEL 


The  iniegnl  (DJ)  is  solved  anmericslly  widi  die  parabolic  role  [Koypers  ad  Tinman,  1963] 
(baaed  on  pair-wiae  q)plkation  of  Sinqison's  mk),  505  intervals  in  first  iteration,  500  increase  for 
eadi  sobeequent  iteratian.  The  integration  limits  sr  a>i  -  (%  -t-  lO-'O  and  oa^  «  1000.  Increasing  die 
vppa  limit  beyond  oo^  ~  1000  did  not  diange  the  value  of  the  integral  significantly.  The  integral 
is  cakulatrd  to  aprecision  of  10-^. 
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IMnJEMENTATIONW  THE  CALCULATION  OTW(xa)  E.1 


PATTON  OF  THE  CALCULATION  OF  W(xj.) 


The  wind  fidd  W(x^)  is  calcuUned  in  a  rectangular,  equidistuttly  qwced  grid  of  Nx  hodzcmtal 
grid  points  and  Nz  vertical  grid  points.  AUhoo^  die  actual  calculation  starts  at  die  water  surface, 
it  is  convoiieiit  to  first  look  at  aregular  grid  p<»nt  wdl  above  the  water  surface.  The  grid  qwcings 
are  denoted  2dx  and  2Az  (see  figure). 


W(x*,zi) 


U(Xj-Ax,Z|-Az) 


U(xrt-Ax,Zi-Az) 


W(xi,Z(-2Ae) 


The  value  of  W  at  a  specific  grid  point  (xi,zj)  is  then  given  1^: 
W(Xi,2i)  =  W(Xi,Zi-2Az)  +  Aw 


(E.1) 


where 


_  U(Xi+Ax,  z^-Az)  -  U(Xi-Ax,z,-Az) 

AW - ^ -  (Ei, 

L^  us  now  turn  our  attention  to  the  first  grid  point  (x^,!^  that  is  at  a  height  wAr  ^  2Az  above  die 
water  surface  T|(Xi).  Close  to  die  wata  surface,  the  gradimts  in  U  and  W  become  huge  and  a  more 
dense  grid  is  needed  to  correcdy  calculate  AW.  Therefore,  «At  »  divided  in  10  equidistant  inter¬ 
vals  widi  qiacing  Ajz.  The  horizontal  grid  is  eiqianding  widi  heigjit: 
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IMPLEMENTATION  OF  THE  CALCULATION  OP  W(x.z) 


X  »  k  A^z 


1  i  k  10 


adieie  k  «  1  in  die  z-interval  ckMest  to  die  surface.  For  eadi  of  the  10  intervals  AW  is  calculated 
and  added  to  die  surface  value  W(Zi,Ti(X|))  to  obtain 

A  difScnlty  arises  when  one  of  the  grid  points  U(xiiAx,  Zi-Az)  used  to  calnilate  AW  is  below  the 
water  surface  (see  figure  above).  In  that  case  the  grid  point  is  neglected  and  AW  is  calculated 
from  the  remaining  grid  point 


AW  = 


±  U(Xi±Ax,  z^-Az)  TU(Xi,Zi-Az) 
Ax 


The  large  number  of  Nz  that  is  needed  for  an  accurate  calculation  of  W(x;z)  presents  anodier 
numerical  problem  since  die  database  of  die  wind  field  W(x«z)  grows  very  large  when  both  Nx 
and  Nz  are  large.  To  keqp  its  size  within  accqitable  limits.  W(x.z)  is  only  calculated  from  the 
wave  trough  to  1  mtr  above  the  wave  crest  Even  then,  it  is  necessary  to  store  only  every  n*^  z-grid 
point  vdiich  results  in  a  z-grid  qiacing  of  typically  1-S  cm  in  die  database  (dqiending  on  die 
wave  heighO-iO 

The  database  may  be  used  to  infer  W(x.z)  at  any  location  (x.z).  To  diis  end,  the  database  is 
searched  frv  the  4  entries  near  (x,z),  see  figure: 


W(x„Zh.i) 


W(x*,.„Zi^,) 


W(x„Zi) 


W(X4^„Zj) 


Note  W(x.z)  descad*  oa  te  kqr  tapot  wUHe.  Thenfcn.  SaaiMfM  of  WOi,^  m  only  valid  te  one  vecifie  vatae  of  Uio- 
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IMnJEMENTATlONOT  THE  CALCULATION  WW(M)  E.3 


Wben  one  of  die  gfid  points  (X|^  or  is  in  die  w«er,  W(x^)  is  cakalsled  from  the  watei 

surfime  (x,T)(x))  m  die  maimer  described  above,  frt  die  odier  case,  first  W(x,Zi)  is  calcolated  by  2- 
point  interpolation  in  die  x-dfrectkn.  Sobseqnently,  W(x,z)  is  calculated  from  W(x,Zi)  by  die 
cmitinuity  equadon.  This  tnediod  yields  accurate  values  for  W(x,z)  a^ieii  die  x-giid  is  densely 
qMced. 

Alternatively,  but  less  accurately,  a  4-point  interpolatian  may  be  used  to  determiiie  W(x,Zj), 
W(x,z^.i)  and  W(x,z),  neglecting  any  missing  grid  points.  Thus,  W(x,Zj)  *  W(xj,Zi)  wbea  (x^.,,Z|) 
is  in  die  water.  However,  in  the  model  code  diis  mediod  is  prefened,  because  it  saves  ccMtqiutation 
time  and  yields  nearly  identical  droplet  trajectories  (see  (4.2). 
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